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Jujuy - Argentina 


Imagen de satélite Landsat TM, en 3D, 
resolución espacial 30 metros, en 
composición Falso Color (5,4,3 / RGB) de la 
Laguna de Pozuelos, Jujuy. Está situada 
sobre una planicie a 3.600 m de altura y 
constituye uno de los humedales más 
importantes a nivel internacional, incluida en 
la lista de la convención de Ramsar. Esta es 
una laguna permanente pero de salinidad 
fluctuante. 
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RESUMEN 


La provincia de Tucumán es la principal productora de 
caña de azúcar de Argentina. El desarrollo de la 
industria azucarera ha hecho que sea la provincia 
mas densamente poblada del país. La quema pre y 
postcosecha de caña es una práctica común, con 
potenciales efectos negativos sobre el ambiente, la 
infraestructura, y la salud de mas de un millón de 
habitantes. En el año 2005 se aprobó una legislación 
prohibiendo la quema de cañaverales; pero hasta el 
momento no se han desarrollado métodos de 
monitoreo para evaluar y efectivizar la legislación. En 
este trabajo presentamos un análisis de focos de 
fuegos en cultivo de caña de azúcar entre los años 


2004-2008, basado en la detección de anomalías 


térmicas por la iniciativa FIRMS y la confección de un 
mapa anual del área cultivada con caña a partir de 
imágenes satelitales TERRA/MODIS. El número de 
focos de fuego descendió del 2004 al 2005, pero 
aumentó sostenidamente en los años subsiguientes, 
alcanzando el máximo valor de la serie (801 focos) en 
2008; lo que implica que la legislación vigente o su 
implementación son ineficientes. Los departamentos 
Cruz Alta, Leales y Simoca concentran el 60% de los 
focos detectados. Agosto y Septiembre concentran la 
mayor densidad de fuego, en coincidencia con las 
menores precipitaciones y humedad atmosférica, lo 
que podría potenciar la contaminación atmosférica y 
el impacto sobre enfermedades respiratorias. 
Finalmente describimos un modelo de monitoreo a 
partir de la metodología propuesta y el sistema de 
modelado grafico STELLA. 


Palabras clave: Imágenes Terra! Modis, 
Contaminación ambiental, Caña de azúcar, 
Tucumán, Fuego 


ABSTRACT 


The Tucumán province is the main sugar cane 
production area in Argentina. The sugar industry has 
made Tucumán the most densely populated province 


in the country. The pre-and post-harvest burning of 
sugar cane is a common practice, with potentially 


negative effects on the environment, infrastructure, 
and health of over one million inhabitants. In 2005, a 
law banning sugar cane burning was passed, but to 
the present there are no monitoring methods to 
evaluate and enforce the legislation. In this work we 


present an analysis of fires hot-spots in sugar cane 


cultivated land between 2004 and 2008. We used 
thermal anomalies by the FIRMS initiative to map fire 


hot-spots, and TERRA / MODIS images to mapps of 


sugar cane cultivated area. The number of fires hot — 


spots fell from 2004 to 2005, but rose steadily in 
subsequent years, reaching the maximum value of 
the series (801) in 2008. This pattern implies that 
existing legislation or its implementation is not 


efficient. Sixty percent of the fires hot-spots were 


detected in departments Cruz Alta, Leales and 
Simoca. 
August and September concentrate the greatest 


number of fire hot-spots, coinciding with the lowest 


rainfall and humidity, which could increase air 
pollution and the effects on respiratory diseases. 
Finally we describe a model of monitoring based in the 
remote sensing and the graphic modelling system 
(STELLA). 


Keywords: Terra! Modis, Environmental pollution, 
Sugar cane, Tucumán, Fires 


INTRODUCCIÓN 


En el Noroeste Argentino, el cultivo de la caña de 
azúcar dio origen a un sistema agroindustrial de gran 
importancia económica y social en las provincias de 
Tucumán, Salta, Jujuy. Particularmente en Tucumán 
la producción cañera fue el motor del desarrollo 
económico. Dos décadas antes de la llegada del 
ferrocarril (1876) la producción de azúcar ya era 
importante, y partir de este año, la superficie de caña 
de azúcar cultivada comenzó a crecer acele- 
radamente al vincularse eficientemente con la 
demanda nacional; y entre 1875 y 1915 la superficie 
aumentó de 2.045 a 106.800 has (Campi, 2000). 
Actualmente existe una superficie cultivada de 
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aproximadamente 220.000 has, según datos de la 
Estación Experimental Obispo Colombres (EEAOC). 
La producción cañera genera múltiples beneficios 
económicos para la provincia, pero conlleva algunos 
problemas ambientales asociados, incluyendo la 
quema de biomasa. La quema de caña previo a la 
cosecha, es una actividad común en muchas áreas 
cañeras del mundo (Hernández, 1995), cuyo principal 
objetivo es aumentar el rendimiento de azúcar en 
fábrica al eliminar residuos a bajo costo y acortar el 
tiempo entre corte y molienda. (Domínguez 
Domínguez, E. et al., 2001; Hernández, 1995; 
Pohlan, H.A.J. et al., 2005). También se lo utiliza con 
la finalidad de disminuir el tiempo entre corte y 
molienda en el ingenio, sin perdida en el contenido de 
sacarosa (Hernández, 1995). Adicionalmente, en 
casos en que la caña se ha cosechado sin quemar, 
también se practica la quema del barbecho para 
favorecer la rápida incorporación de materia orgánica 
al suelo y para evitar que quemas accidentales del 
barbecho en los meses posteriores afecten el rebrote 
de la cepa de caña correspondiente a la siguiente 
campaña agrícola. Por otra parte la quema de caña 
tiene consecuencias negativas sobre el medio 
ambiente: disminuye el contenido de nitrógeno, la 
porosidad y humedad del suelo; aumentando su 
compactación y susceptibilidad a la erosión 


(Dourado—Neto, D. etal., 1999; Oliveira, J.C.M. de, et 


al. 2000; Tominaga, T.T. et al., 2002). Además la 
quema propicia la pérdida de nitrógeno del suelo 
(Domínguez Domínguez et al., 2001). Pero 
posiblemente la consecuencia ambiental más 
importante sea el efecto sobre la contaminación 
atmosférica. La quema libera CO2, CO y O3 en 
cantidades significativas (Kirchoff, V.W.J.H. et al., 
1991) y abundantes partículas sólidas que tienen 
consecuencias para la salud humana y la visibilidad, 
generando accidentes y afectando el turismo. 
Adicionalmente, el efecto directo del fuego (e.g. 
quema de generadores de luz y otra infraestructura) 
tiene fuertes impactos económicos. En respuesta a la 
importante percepción social de estos problemas, se 
sancionó en 2005 una modificación del artículo 38 de 
la ley provincial n* 6253 que prohíbe la quema de 
cañaverales en todo el territorio de la provincia. Sin 
embargo, hasta el presente no se cuentan con 
mecanismos de monitoreo sistematizados que 
permitan evaluar el efecto de la ley y promover su 
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aplicación. 

Los sensores remotos basados en satélites proveen 
herramientas altamente eficaces y de bajo costo para 
este propósito. La teledetección ha sido empleada 
para el monitoreo de fuegos desde la década de los 


90s basándose en el satélite NOAA-AVHRR y más 


recientemente en el satélite TERRA/MODIS, en 
particular su programa de detección de anomalías 
térmicas que ha comenzado a ser utilizado 
regularmente para detección de fuegos a grandes 
escalas espaciales (Bowman, D. M. J.S. et al., 2009) 
a la continental (Di Bella, C.M. et al., 2006; Chuvieco, 
E.etal., 2008). 

Las imágenes TERRA MODIS han sido usadas 
recientemente para la estimación de la superficie 
cultivada con caña de azúcar en Tucumán (Soria, F.J. 
etal., 2000; Soria, F.J. et al., 2007) y diversos cultivos 
del noroeste argentino (Volante, J. et al. , 2004). Sin 
embargo, y a pesar de la importancia que tiene la 
quema del cultivo de caña de azúcar para la provincia 
de Tucumán, hasta el momento no se dispone de 
estimaciones cuantitativas de este proceso ni 
indicadores de cuales son sus patrones temporales y 
espaciales. En este trabajo presentamos un 
seguimiento de los focos de fuego de la provincia de 
Tucumán asociados al cultivo de caña de azúcar 


durante el periodo 2004-2008, basado en la iniciativa 


"Fire Iniciative for Resources Managemente 
Systems” (FIRMS, 2009). 


MATERIALES Y MÉTODOS 
ÁREADE ESTUDIO 


La provincia de Tucumán se encuentra ubicada en el 
Noroeste de la Argentina, entre los S 26 y S 28 ? latitud 
y entre los O 64.5" y O 66" de longitud. Los principales 
cultivos se localizan en las zonas de pedemonte 
(caña de azúcar, limón, hortalizas) y la llanura del 
Este de la provincia (soja, maíz, trigo). El clima es de 
tipo monzónico con la estación seca durante el 
invierno y una estación de lluvia durante el verano. La 
caña de azúcar se cultiva por encima de los 800 mm 
de precipitaciones. El área de cultivo encuentra un 
límite en el este asociado con las lluvias y al oeste con 
la sierra de Aconquija que presenta suelos con altas 
pendientes (Figura 1) 
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Figura 1. Área de estudio. Distribución de | cultivo de caña de azúcar en los distintos departamentos de la provincia de 


Tucumán. 


CARTOGRAFÍA DE LA SUPERFICIE DEL CULTIVO 
DE CAÑA DE AZÚCAR PERIODO 2004-2008 


Para restringir el análisis a fuegos asociados a 
cultivos de caña de azúcar realizamos un mapa de 
este cultivo en la provincia de Tucumán para cada uno 
de los años estudiados. Para confeccionar el mapa se 
empleó el producto de índices de vegetación del 
satélite Terra/Modis (MOLT/MOD13Q1.005/VI 
_T6dias_250). Este producto es una composición de 
observaciones diarias para un periodo de 16 días, 
tienen una resolución espacial de 250 m e incluyen 
cuatro bandas (Azul, Rojo, Infra Rojo e Infra Rojo 
Medio) además de dos índices de vegetación: el 
Índice de Vegetación Normalizado (NDVI por sus 
siglas en ingles) y el Índice de vegetación realzado 
(EVI por sus siglas en ingles). En este trabajo 
empleamos las imágenes que se componen entre los 
días 7 y 23 de Abril por considerarse que es la fecha 
óptima para poder realizar una buena detección del 
cultivo al corresponder a un momento donde la caña 
esta bien desarrollada y aún no se inicio la cosecha. 
Las imágenes fueron acondicionadas para este 
trabajo llevándolas al sistema de coordenada oficial 


de Argentina Gauss-Krugger (Sistema POSGAR94) 


Faja 3. No se realizaron correcciones geográficas o 
atmosféricas adicionales ya que el producto es 
proporcionado por la NASA con las calibraciones 
correspondientes. De la escena completa se realizó 
un recorte para restringir los datos al área de la 
provincia de Tucumán. Para el procedimiento de 
clasificación se utilizó la información espectral de las 
cuatro bandas a las cuales se le agregó la altitud 
como una banda más; derivada de un modelo de 
elevación digital DEM, htpp://SRTM.CSI. 
CGIAR..ORG. 

La inclusión del DEM es de utilidad para diferenciar el 
cultivo de caña ya que este cultivo se restringe a las 
zona de llanura y respeta un limite altitudinal claro, por 
lo cual la altitud permite diferenciar tipos de 
vegetación con firmas espectrales similares a la caña 
pero ubicados en altitudes diferentes (e.g. pasturas). 
Se realizó una clasificación supervisada definiendo 
tres clases (caña, suelo desnudo y bosque). Las 
áreas de entrenamiento fueron distribuidas por toda 
el área de estudio y se utilizaron al menos 70 sitios 
para cada clase. Estas clases de entrenamiento se 
definieron a base de conocimiento del terreno. El 
clasificador empleado fue el método Maximun 
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likelinood con una probabilidad del 0.9 en base al 
programa ENVI 4.55. Debe señalarse que este 
procedimiento no pretende hacer un mapa detallado 
de tipos de vegetación si no disponer de una 
estimación particular de la zona cultivada de caña de 
azúcar. 


DETECCIÓN DE FOCOS DE FUEGO 


Los mapas de fuegos activos de la provincia de 
Tucumán fueron obtenidos a partir de la detección de 
focos de fuegos activos provistos por FIRMS (2009). 
Este sistema integra información proveniente de 
sensores remotos con Sistemas de Información 
Geográfica (GIS) con el fin de proveer mapas 
globales de ubicación de anomalías térmicas (focos 
de fuegos/puntos de calor) detectados a partir de 
imágenes satelitales MODIS (Moderate Resolution 
Imaging Spectroradiometer —ftp://e4ftl01u.ecs.nasa 
.gov/MOLT/MOD13Q1.005). La ubicación de los 
focos de fuegos es procesada por el Sistema de 
Respuesta Rápida de MODIS/NASA utilizando el 
producto estándar de MODIS “MOD14 Fire and 
Thermal Anomalies Product”. El algoritmo de 
detección de fuegos/puntos de calor, se basa en la 
información contenida en las bandas de infrarrojo 
medio e infrarrojo térmico de MODIS. En la mayoría 
de los casos, la anomalía térmica detectada de esta 
manera se debe a un fuego pero a otras fuentes de 
calor como actividad volcánica, chimenea de venteo 
de pozos de gas, etc. En condiciones ideales 
(observación nadir, superficie homogénea, escena 
libre de nubes y libre de humo), el algoritmo puede 
llegar a detectar fuegos de 50 m2; pero en 
condiciones mas generales puede suponerse una 
detección confiable de áreas quemadas de hasta 100 
m2. Cada foco de fuego detectado se asigna a un 
pixel de 1 km x 1 km que en el algoritmo representa la 
ocurrencia de uno o más fuegos/puntos de calor. Si 
bien el “evento de fuego” se asigna a todo el píxel, 
típicamente los fuegos detectados sólo afectan una 
fracción del mismo, por lo que el método no es preciso 
para cuantificar área en una determinada fecha 
(asignar el fuego a toda la superficie del píxel 
implicaría una sobreestimación sistemática del área 


de fuegos). Sin embargo es una medida no-sesgada 


y por lo tanto confiable de variabilidad temporal y 
permite hacer comparaciones confiables de cantidad 
de fuegos entre distintos años. 

Los fuegos detectados por el sistema FIRMS, 
presentan una estimación de la confiabilidad y 
registra el día en el cual fue detectado. En este trabajo 
consideramos los focos con una confianza superior al 
50%. Por último, utilizando el mapa de área cultivada 
de caña se seleccionan los focos de fuegos que 
ocurren en sitios con caña de azúcar al inicio de la 
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cosecha (Abril) con lo cual se restringen los focos de 
fuegos asociados a este cultivo eliminando otras 
situaciones como quema de bosque para desmonte. 


RESULTADOS 


Según nuestras estimaciones la superficie de caña de 
azúcar se redujo en aproximadamente 40000 ha 
entre 2004 y 2005 (20% del total) para volver a 
ascender y estabilizarse en los años subsiguientes en 
poco mas de 200000 ha (Figura 2). Nuestras 
estimaciones de superficie de caña de azúcar 
siempre tuvieron una discrepancia menor del 10% 
con respecto a las estimaciones mas detalladas 
realizadas por la EEAOC, por lo cual se puede 
considerar que nuestras estimaciones de áreas 
cultivadas son aceptables para los objetivos de este 
trabajo. 

La distribución temporal de fuegos también mostró 
una reducción marcada en el año 2005, en forma 
coincidente con la aprobación del decreto N* 795/3, 
que modifica el articulo 38 de la ley N*6253, y prohíbe 
la quema de caña de azúcar como método auxiliar 
para la cosecha; y con la reducción de la superficie de 
caña de azúcar. Sin embargo, entre el 2006 y 2008 el 
número de fuegos aumenta sostenidamente (en 
promedio cercano a un 30% de incremento anual), 
con lo que el número de fuegos detectados en 2008 
(801) duplica a los detectados en 2005. Alo largo de la 
serie el número de focos de fuego acompaña 
aproximadamente la curva de superficie de caña pero 
es mucho más abrupta. 

La dinámica de focos de fuego durante los años del 
estudio (Figura 3) muestra que lo primeros fuegos 
asociados al cultivo de caña de azúcar ocurren entre 
los meses de Mayo y Junio, lo cual coincide con el 
comienzo de la zafra en la provincia de Tucumán; y 
alcanzan un pico en los meses de agosto y 
septiembre. 

Durante el período de estudio el 95% de los fuegos 
ocurrieron entre Julio y Octubre, siendo Septiembre el 
de mayor porcentaje (29%). La ocurrencia de fuegos 
coincide precisamente con la baja ocurrencia de 
precipitaciones (Figura 3). 

En cuanto la distribución espacial de los focos de 
fuegos (Figura 4) en los diferentes departamentos de 
la provincia de Tucumán durante los cinco años 
estudiados, se puede observar que los 
departamentos de mayor superficie de caña de 
azúcar, Cruz Alta, Leales y Simoca, son los que 
presentan los mayores porcentajes reuniendo entre 
los tres el 60% del total de fuegos monitoreados en 
estos 5 años (Figura 5). Los restantes Departamentos 
presentan valores inferiores al 8% del total de fuegos 
en el período 2004 — 2008. 
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Figura 2. Superficie estimada de caña de azúcar y focos de fuego detectados en el área cañera de la provincia de 
Tucumán entre 2004 y 2008. 
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Figura 3. Distribución mensual de precipitaciones (Fuente: SAGPYA) y ocurrencia de fuegos en el área cañera 
de la provincia de Tucumán entre 2004 y 2007. 
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Figura 4. Focos de fuego sobre las áreas cultivadas con caña de azúcar en la provincia de Tucumán, para el 
periodo 2004-2008. 
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Figura 5. Departamentos de Tucumán con mayor cantidad de focos de fuegos durante el período 2004-2008 
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DISCUSIÓN Y PROYECCIONES A FUTURO 


Este trabajo presenta la primera cuantificación de 
patrones temporales de fuego a lo largo de un período 
de varios años en el área cañera de la provincia de 
Tucumán. El sistema de monitoreo puede jugar un rol 
central en evaluar la eficacia de la ley y sus 
mecanismos de implementación. De hecho, la 
difusión de resultados preliminares de este trabajo ha 
tenido un fuerte impacto mediático, motivando 
importantes operativos de control durante 2009 (La 
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Gaceta, 2009). Dada las fuertes implicancias 
ambientales, el monitoreo de fuegos puede servir 
como insumo central para un modelo integrado de 
evaluación permanente de las consecuencias 
ecológicas de la quema de cañaverales, que incluya 
los costos por emisión de gases de invernadero, los 
efectos en la salud (principalmente enfermedades 
respiratorias), la pérdida de fertilidad, y las 
consecuencias para la infraestructura y accidentes 
derivados de falta de visibilidad en rutas (Figura 6). 
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Figura 6. Esquematización de un modelo integrado de monitoreo y evaluación permanente de los efectos ambientales de 
los fuegos de cañaverales, en lenguaje de modelado gráfico STELLA. Recuadros simbolizan variables acumulables como 
partículas en la atmósfera o nitrógeno en el suelo; líneas dobles con círculos anexos simbolizan flujos de materia; 
circulos simbolizan variables de control; líneas delgadas simbolizan flujos de información. 
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En el esquema se representan los principales flujos 
asociados a la quema de cañaverales. Los fuegos 
influyen en las emisiones de COZ2 a la atmósfera, en la 
acumulación de partículas sólidas en la atmósfera, en 
las pérdidas de Nitrógeno del suelo y en la facilitación 
de accidentes y destrucción de infraestructura. Todos 
estos procesos generan distintos costos que suman a 
un costo total del proceso. 

En este trabajo la superficie cultivada con de caña fue 
estimada mediante clasificación de imágenes 
satelitales cuyo mapas resultantes no disponen de 
una verificación de campo pero coincide con 
estimaciones independientes. Los mapas obtenidos 
si bien pueden presentar errores son lo 
suficientemente precisos para nuestro propósito de 
detectar fuegos asociados con la producción de caña. 
Por otra parte debe indicarse que no hay en la zona 
otros cultivos que espectralmente se puedan 
confundir con caña de azúcar y que sean manejados 
con fuego. Por lo tanto los posibles errores de los 
mapas del área cultivada con caña son menores y no 


ocasionan sobre-estimación de fuegos de caña ya 


que si se incluye como caña erróneamente algún otro 
cultivo o tipo de vegetación estos sitios no están 
sometidos a fuego. Por el contrario los errores de 
omisión del mapa (es decir la no detección de áreas 
cultivadas con caña) puede significar una 
subestimación de los fuegos. Teniendo en cuenta lo 
anterior, junto con la utilización de un umbral de 
confianza del producto de fuegos del 50% y que un 
píxel de 1 km indicado como una anomalía térmica, 
puede incluir más de un fuego activo nosotros 
argumentamos que nuestras estimaciones de 
números de fuego son conservadoras. Con todas 
estas limitaciones, lo mas importante de destacar es 
que el método no es sesgado en el tiempo, y permite 


una cuantificación objetiva de patrones inter-anuales 


y estacionales, que son esenciales para un sistema 
de monitoreo y su uso para políticas ambientales. 
Para mejorar la calibración y precisión del método en 
el caso de la quema de caña de Tucumán algunos 
preguntas que deberían resolverse incluyen: ¿Cual es 
el umbral de la confiabilidad adecuado para detectar 
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fuegos? ¿Cual es la dimensión mínima de un fuego 
en el cultivo de caña de azúcar, para que sea 
detectado por el satélite? ¿Ocurren eventos de 
fuegos que por tener corta duración, no están activos 
en el momento de registro de la imagen y por lo tanto 
no se detectan? ¿Como estimar superficie quemada 
a partir de fuegos registrados en píxeles de 1 Km2? 


CONCLUSIONES 


* El mapeo de fuegos basado en la detección de 
anomalías térmicas por la iniciativa FIRMS, 
combinado con mapas anuales del área cultivada con 
caña a partir de imágenes satelitales TERRA/MODIS, 
permiten el desarrollo de un sistema de monitoreos 
de fuegos en cultivos de caña de azúcar, objetivo, de 
buena resolución temporal y espacial, altamente 
confiable y de bajo costo; que podría integrarse con 
modelos dinámicos que cuantifiquen las distintas 
consecuencias ambientales del proceso, y sus costos 
para la sociedad. 

* Entre 2004 y 2008 se observó una disminución 
inicial de los fuegos en el año 2005, coincidente con la 
sanción de la prohibición provincial del uso de fuegos, 
y con la disminución del área cañera. A partir de esa 
fecha, los fuegos han aumentado sostenidamente, 
evidenciando la ineficacia de la ley o sus mecanismos 
de implementación. 

e La mayor ocurrencia de focos de fuegos es durante 
la estación seca con lo cual se acentúan los 
problemas ocasionados por la contaminación 
atmosférica, tales como efectos en la visibilidad o 
incidencia en enfermedades respiratorias. 
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RESUMEN 


El fuego es una herramienta frecuentemente utilizada 
como parte del manejo agropecuario. Su uso 
indebido puede constituir una práctica riesgosa, 
generando efectos no deseados con daños a la 
producción agropecuaria y el ambiente. En 
teledetección se utilizan generalmente sensores 
térmicos y ópticos para la identificación de puntos de 
calor y la estimación de daño por fuego en la 
vegetación, respectivamente. Los sistemas de 
RADAR pueden presentar ventajas dada su 
capacidad de penetración en un canopeo y la 
posibilidad de obtener imágenes en condiciones de 
nubosidad o con presencia de humo. El RADAR en 
banda L posee mayor capacidad de penetración en la 
vegetación, pero ha sido escasamente estudiado 
para el análisis de eventos de fuego. En este trabajo, 
se consideró un evento de fuego ocurrido en la región 
del delta del río Paraná durante el año 2008. Se 
realizaron comparaciones entre clasificaciones de 
área quemada obtenida con el RADAR ALOS 
PALSAR (banda L) y clasificaciones obtenidas con 
sensores ópticos y térmicos. Se consideró también el 
resultado de la comparación según diferentes tipos 
de vegetación. Se analizó la clasificación con ALOS 
utilizando umbrales de cambio en el nivel de 
retrodispersión. Se observó una asociación entre los 
resultados obtenidos para los sensores de RADAR y 
ópticos mayor al 70 %. Considerando la imagen en su 
conjunto, se estableció un valor umbral óptimo para la 
clasificación con ALOS, alrededor de 5dB. Diferentes 
valores umbrales pueden ser sugeridos para los 
distintos tipos de vegetación considerados. 


Palabras clave: Radar, Landsat, Vegetación, Área 
quemada, Fuego 


ABSTRACT 


Fire is commonly used as part of agricultural 
practices. The inadequate use of this tool can be 
dangerous, generating undesired effects to crop and 
livestock production and to the environment. 
Common remote sensing tools for hotspots detection 


Argentina 


and burnt vegetation damage are based on thermal 
and optical sensors. RADAR systems can take 
advantages to these techniques, since its ability to 
penetrate vegetation canopies and the ability to obtain 
images despite the presence of clouds or smoke 
derived from fires. L band RADAR systems can 
penetrate greater in vegetation systems, but have 
been scarcely used for fire events analysis. In this 
work we analyzed the fire event occurred at the 
Paraná river delta region on 2008. 

Comparisons between burnt area classifications 
obtained from L band ALOS PALSAR RADAR, optical 
and thermal sensors were performed. Also 
comparison results were performed taking into 
account different vegetation units. Comparison 
results were analyzed using different backscatter 
change thresholds for ALOS images. The results 
showed an association between RADAR and optical 
based classifications of over 70 %. For the whole 
image, an optimal threshold value of 5 dB was 
determined. Also different threshold values can be 
suggested for different vegetation units. 


Keywords: Radar, Landsat, Vegetation, Burnt area, 
Fire 


INTRODUCCIÓN 


El fuego es una herramienta comúnmente utilizada 
como parte del manejo agropecuario en distintas 
zonas de Argentina (Kunst et al., 2003). Idealmente 
bajo un esquema de quemas prescriptas oO 
controladas, el objetivo de quemar está puesto en la 
obtención de material palatable para el ganado, la 
reducción de combustibles, la eliminación de 
residuos de explotación forestal y la ampliación de 
nuevas tierras productivas, entre otros usos (Kunst et 
al., 2003). A pesar de los beneficios que se pueden 
obtener del uso del fuego, además de la concepción 
del mismo como un disturbio regulador de ciertos 
ecosistemas, su utilización desmedida puede 
constituir una práctica riesgosa capaz de producir 
efectos no deseados tanto desde el punto de vista 
productivo y antrópico a corto plazo, como ambiental 
a largo plazo. 


1 > Edición Especial Vol. 2 


piciÉmBRE 2010. EVALUACIÓN DEL USO DE IMÁGENES DE RADAR EN BANDA L 


Uno de los aspectos de mayor relevancia que atañen 
al uso del fuego como herramienta es la escala con la 


cual se producen las quemas (local-regional-global), 


la superficie quemada, la intensidad y la frecuencia. 
Estos elementos se corresponden con algunos de los 
atributos del régimen de fuegos de determinada 
región, que a partir de la intervención del hombre ya 
sea por el uso indiscriminado o la supresión de 
fuegos, se ven notablemente alterados, con las 
consecuencias que ello implica en la estructura y el 
funcionamiento de los ecosistemas, las alteraciones 
en la composición del suelo, los efectos producidos 
por la emisión de gases de efecto invernadero hacia 
la atmósfera además de las perdidas a nivel humano 
y económico. 

En términos de prevención, cuando una quema se 
torna incontrolable se convierte en un incendio, y de 
acuerdo a la escala y los recursos en riesgo pueden 
convertirse en un desastre natural. En este contexto, 
existe la necesidad de registrar el desarrollo de las 
quemas y su evolución hasta convertirse en un 
desastre natural. 

Las técnicas de teledetección principalmente 
utilizadas para identificar fuegos activos, como para 
determinar la superficie quemada y estimar el daño 
en la vegetación producido por el fuego, están 
relacionadas con sensores térmicos y ópticos 
respectivamente o bien a partir de una combinación 
entre ellos (Chuvieco et al., 2008; Fraser et al., 2000; 
Huesca etal., 2008). 

La identificación de un fuego activo o punto de calor 
implica por lo general establecer umbrales en la 
temperatura de brillo captada por el sensor en torno a 
los 4 micrones, y en la diferencia entre los 4 y 11 
micrones (Giglio et al., 2003). Actualmente el sensor 
más utilizado para tal fin es el MODIS, cuyo diseño 


proviene de la experiencia del sensor NOAA-AVHRR 


y GOES (Justice et al., 2002). La detección de Áreas 
quemadas por medios ópticos se realiza 
comúnmente a partir de los canales visibles, Infrarrojo 
cercano y el infrarrojo medio de onda corta (SWIR) 
especialmente sensibles a la señal del carbón 
(Pereira 1999). Existen diversas técnicas que 
permiten obtener los perímetros quemados a partir de 
dichos canales que van desde la interpretación visual, 
hasta la aplicación de índices especificamente 
diseñados, como el NBR, dNBR BAIM, NDVI, entre 
otros. La estimación de la superficie y el nivel de daño 
se realizan tanto a partir de imágenes de moderada 
resolución (MODIS) o alta resolución (LANDSAT, 
CBERS) para cartografías de mayor detalle. 

Las imágenes de RADAR pueden proporcionar 
información adicional para la estimación de daño por 
fuego, debido tanto a su capacidad de extraer 
imágenes con presencia de nubosidad o humo, como 
a la capacidad de penetración de la señal emitida 
dentro del canopeo. Esta última característica 
permite obtener información sobre la estructura de la 
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vegetación (Le Toan et al., 2004), pudiendo llegar a 
brindar información sobre la biomasa vegetal 
presente. En Argentina, se ha evaluado el uso de 
RADAR en banda C (5.3 GHz) en el delta del río 
Paraná para la a estimación de áreas quemadas de 
juncales y pajonales (Karszenbaum etal., 2003). 

El uso de RADAR en banda L (1.27 GHz ) presenta 
mayor capacidad de penetración en canopeos 
vegetales, mostrando menores inconvenientes de 
saturación de su señal con altos niveles de biomasa 
vegetal. 

Actualmente existe un único satélite en órbita que 
dispone de un RADAR en banda L(ALOS PALSAR) y 
su utilización ha sido escasamente evaluada en 
Argentina. 

Se han realizado algunos trabajos para la 
caracterización de áreas quemadas con imágenes de 
RADAR mediante el análisis del valor de 


retrodispersión entre imágenes pre y post-fuego 


(Honikel y Wagemuller, 2000; Karzenbaum et al., 
2003; Shimada et al., 2000). Sin embargo, es difícil 
encontrar en la bibliografía valores umbrales de 
cambio en los niveles de retrodispersión para 
estimaciones de área quemada, particularmente para 
banda L, lo que facilitaría la realización de 
clasificaciones a partir de imágenes de RADAR. La 
capacidad de obtener información polarimétrica 
(característica incorporada en los satélites más 
recientes: ALOS PALSAR, RADARSAT 2, TerraSAR), 
permitiría obtener mayor información sobre 
características de la vegetación. En particular se 
considera que la polarización donde se envía la señal 
en dirección horizontal y se recibe su rebote en 
dirección vertical (HV), está muy relacionada con 
características de la vegetación. 

Este trabajo se centra en el evento ocurrido en la 
región del delta del río Paraná durante el año 2008. 
Allí se produjeron incendios que afectaron una región 
muy importante del delta e incluso, debido al humo 
producido, se sintieron efectos dañinos en centros 
poblacionales cercanos. Si bien es difícil conocer la 
causa de eventos de fuegos, estos coincidieron con 
un período de fuerte sequía observado en la región 
durante 2008, que fue antecedido por un período muy 
húmedo con importantes inundaciones durante 2007. 
Esta situación pudo haber generado grandes niveles 
de biomasa vegetal combustible, y condiciones 
predisponentes tanto para el inicio de los fuegos 
como para su propagación posterior. El mayor 
número de puntos de calor detectados se produjo en 
el mes de abril de 2008, siendo importante también la 
aparición de puntos dentro del período de marzo a 
octubre de 2008 (Liota, 2008). 

El objetivo de este trabajo es la evaluación de 
técnicas para la detección de daño por fuego en 
vegetación obtenidas a partir de imágenes de 
RADAR en banda L, mediante la comparación con el 
resultado de técnicas comúnmente utilizadas a través 
de sensores ópticos y térmicos. 
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METODOLOGÍA 
ÁREADE ESTUDIO 


El delta del río Paraná se forma a partir de la 
acumulación de sedimentos en una llanura de 
inundación formando una gran cantidad de islas. 
Presenta una gran heterogeneidad de unidades de 
paisaje, pudiéndose encontrar, particularmente en la 
región de análisis de este trabajo, desde grandes 
zonas con relieve elevado y presencia de praderas 
graminiformes (isla de lbicuy), hasta zonas 
inundadas en forma permanente y no permanente 
(Malvárez, 1997). Una de las situaciones más 
frecuentes es la presencia de islas conformadas por 
albardones perimetrales (con presencia de bosques 
naturales e implantados) que rodean áreas 
deprimidas (presencia de praderas de herbáceas, 
pajonales o juncales). 

Se señalan cinco tipos de vegetación principales en la 
zona del delta del Paraná (Salvia etal., 2009): 


- Bosques de sauce (Salix sp.) y álamo (Populus sp.), 
espinillo (Acacia caven) y algarrobo (Prosopis sp.) 
- Pajonales de cortadera (Scirpus giganteus) y 
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Panicum pronitis. 
- Juncales de Schoenoplectus californicus. 


- Praderas de herbáceas acuáticas graminoides y 
latifoliadas 
- Pastizales de gramíneas. 


La forestación se ha convertido en una de las 
actividades económicas principales del delta durante 
las ultimas décadas, a partir de la reconversión desde 
actividades de fruticultura, horticultura y producción 
de mimbre (Borodowski, 2006). La ganadería 
también es una actividad de importancia, que 
adquiere relevancia durante los últimos años a partir 
de la expansión de la frontera agrícola y el 
desplazamiento de la ganadería hacia zonas 
marginales como el delta. 


DETECCIÓN DE ÁREAS QUEMADAS 


Se realizó una comparación entre estimaciones de 
área dañada por fuego en una zona del delta del río 
Paraná obtenida a partir de imágenes de RADAR y 
estimaciones realizadas a partir de sensores ópticos 


(Fig. 1). 
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Figura 1. Área de estudio. Rojo: Perímetro de imagen LANDSAT 5 TM (PATH: 226; ROW: 83); Azul: Perímetro 
de imagen ALOS PALSAR (PATH: 90; FRAME: 6500). Area rayada: Región analizada en este trabajo. 
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IMÁGENES UTILIZADAS 


Se utilizaron dos imágenes del RADAR ALOS 
PALSAR modo FBD (Polarización HH y HV) nivel de 
procesamiento 1.0 correspondientes a la zona del 
delta del río Paraná captadas los días 23/07/2007 
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(preevento) y 25/07/2008 (post evento) (Fig. 1 y 2). 
Como imagen óptica se utilizó una imagen LANDSAT 
5 TM captada el día 06/10/2008 (post-evento) (Fig. 1 
y 3). Se utilizaron también imágenes MODIS 
captadas entre enero y julio de 2008 para la detección 
de puntos de calor. 





Figura 2. a) Imagen ALOS PALSAR FBD del 23/07/2007 (pre evento). Canales: Rojo: HH; Verde: HV; Azul: HH. 
b) Imagen ALOS PALSAR FBD del 25/07/2008 (post-evento). Canales: Rojo: HH; Verde: HV; Azul: HH. 


PROCESAMIENTO IMÁGENES ALOS 


Se utilizó el software RO! pac (Rosen, Hensley, 
Peltzer, 4 Simons, 2004) para el procesamiento de 
las imágenes ALOS PALSAR nivel 1.0 (datos crudos), 
con lo que se obtuvieron imágenes en formato single 
look complex (información de amplitud y fase). Se 
realizó una compresión (multilook) en azimut con 
factor 5 para obtener un tamaño de píxel que 
represente una superficie cuadrada y se estimó la 
intensidad de retrodispersión. Las imágenes fueron 


geo-—referenciadas y corregidas geométricamente 


utilizando como base de comparación a imágenes 
LANDSAT ETM ortorectificadas provistas por el 
programa Global Land Cover Facility 
(http: //glcfapp.umiacs.umd .edu:8080/esdi/index.jsp) 
de la Universidad de Maryland (EEUU). El tamaño de 
píxel finalmente obtenido fue de 18.8 m. Se realizó un 
filtrado con el método FROST utilizando un tamaño 
de ventana de 5 píxeles para reducir el ruido de la 


imagen (speckle). 

Para considerar los cambios en la vegetación entre la 
imagen previa y la posterior al evento de fuego, se 
estimó la diferencia de intensidad de retrodispersión 
entre el canal de polarización HV de las dos imágenes 
de RADAR, dado que se considera que esta 
polarización es mejor indicadora de niveles de 
vegetación. 

Se realizó una clasificación no supervisada con el 
método ISODATA (Jensen, 2000). Se consideraron 3 
clases que representarían tres estados de vegetación 
en relación con la intensidad anterior al evento (2007) 
y posterior al evento (2008): 1) diferencia positiva, 2) 
diferencia negativa y 3) sin cambios. Se consideró 
como área quemada a la clase con diferencias 
positivas, es decir, que presentó mayor nivel de 
intensidad pre-evento (2007) que post-evento 
(2008). Posteriormente se realizó un filtro de mediana 
con un tamaño de ventana de 9 píxeles, para reducir 
la presencia de píxeles aislados. 
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Figura 3. Imagen LANDSAT TM del 06/10/2008 (post evento). Canales: Rojo: Banda 7; Verde: Banda 4; Azul: Banda 3) 


PROCESAMIENTO IMÁGENES LANDSAT 


Se realizaron correcciones geométricas de las 
imágenes utilizando como base de comparación a 
imágenes LANDSAT ETM ortorectificadas provistas 
por el programa Global Land Cover Facility 
(http: //glcfapp.umiacs.umd.edu:8080/esdi/index.jsp) 
de la Universidad de Maryland (EEUU). Los valores 
digitales (digital number) de las bandas 3,4 y 7 de 
LANDSAT 5 TM fueron transformados a irradiancia y 
posteriormente a reflectancia, siguiendo la 
metodología descripta por Chander y Markham 
(2003). A posterior se calculó el índice de área 
quemada NBR (Normalized Burned Ratio) a partir de 
las bandas 4 y Y de acuerdo a la metodología 
desarrollada por Key y Benson (2005) según la 
fórmula: (NBR = R4 — R7 / R4 + R7). La detección y 
posterior delimitación de la superficie quemada se 
realizó de acuerdo a los umbrales sugeridos por Key 
á Benson (2005) y adaptados bajo un ejercicio de 
prueba y error. La cicatriz finalmente obtenida se 
reclasificó a dos únicas clases, quemado y no 
quemado. 


PROCESAMIENTO IMÁGENES MODIS 


Se utilizaron los puntos de calor (Hotspots) MODIS 
provistos por CONAE procesados con el algoritmo 
MODO014 v3.2 (Justice et al., 2006). Se acumularon 
las detecciones correspondientes al período 


comprendido entre las imágenes ALOS. Luego se 
generó un buffer cuadrado de 1 Km de lado, ubicado 
en la coordenada central del punto de calor 
correspondiente al tamaño de píxel MODIS. Se 
generó una imagen binaria entre área cubierta y no 
cubierta por puntos de calor. 


COMPARACIÓN ENTRE CLASIFICACIONES 


Se utilizó una máscara para filtrar áreas agrícolas y 
urbanas y cuerpos de agua, y se consideró una 
intersección de las áreas cubiertas por las imágenes 
utilizadas. El tamaño de píxel de los distintos tipo de 
imagen fue resampleado a 18.8 m (tamaño de píxel 
obtenido para las imágenes ALOS). Como resultado 
de la comparación se establecieron cuatro 
categorías, según el criterio presentado en el cuadro 
1. Se realizaron comparaciones entre imágenes 
ALOS y LANDSAT, entre imágenes LANDSAT y 
Hotspots, y entre imágenesALOS y Hotspots. 

A partir de una clasificación portipos de vegetación en 
la zona del delta del río Paraná (Salvia et al., 2009). 
Se analizaron los resultados las clasificaciones con 
ALOS y LANDSAT para diferentes tipos de 
vegetación. 

En esta clasificación se diferencian cinco tipos: 
juncales, cortadera, herbáceas acuáticas, pastizales 
y bosques nativos y plantaciones forestales. 
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Cuadro 1: Identificación de resultados de la clasificación a partir de la combinación de las distintas clasificaciones. 


Clasificación imagen 1 
Quemado 
Quemado 

No quemado 


NO cuemaco NO 





DETECCIÓN DE UMBRALES ALOS PALSAR 


Se establecieron posibles umbrales de disminución 
de la intensidad de retrodispersión (HW 2007- 


HV2008) para la clasificación con imágenes ALOS. 
Se consideraron valores umbrales entre 2 y 8 dB. Se 
realizaron comparaciones entre clasificaciones con 
ALOS y LANDSAT utilizando los diferentes umbrales 
definidos. A partir de su análisis se determinó un valor 
umbral óptimo de ALOS para la detección de área 
dañada por fuego. Se analizó también la variación de 
estos umbrales según el tipo de vegetación. 


RESULTADOS Y DISCUSIÓN 


Se observa una diferencia importante entre la imagen 
ALOS PALSAR anterior y la posterior al evento (Fig. 
2). En la imagen de 2008, se observa una gran 
proporción de la superficie del delta con colores 
oscuros (escasa intensidad de señal de RADAR), a 
diferencia de la imagen de 2007 donde se observa 
mayor proporción de colores claros y verdosos 
(mayor importancia relativa de canal HV), lo que 
indicaría una mayor presencia de vegetación en 
aquel momento. 

La diferencia en intensidad de retrodispersión del 
canal HV entre la imagen previa (2007) y la posterior 


al evento (2008 - Fig. 4.a) puede ser un indicador de 


cambios en la vegetación. Las zonas de color oscuro 
(valores altos y positivos de la diferencia de 
intensidad), indicarían disminuciones importantes en 
los niveles de vegetación. Es importante tener en 
cuenta que pueden ocurrir efectos confundidos en la 
señal del RADAR debido a cambios en la humedad 
del suelo y de la vegetación entre ambas imágenes, 
particularmente debido a que en la zona se produjo 
una inundación importante durante 2007 y una fuerte 
sequía durante 2008. Cambios negativos en la 
intensidad indicarían incrementos en el nivel de 
vegetación (colores claros — Fig. 4.a), pero también 
esto podría observarse en zonas que estuvieron 
inundadas durante 2007 (la señal es reflejada) y no lo 


Clasificación Imagen 2 
Quemado 

No quemado 
QCuemaco 

vemado 





ñesultado comparación 
ACIERTO QUEMADO 
OMISON 
LOMBION 
ACIERTO NO QUEMADO 





estén durante 2008 (se capta señal del RADAR). Los 
colores grises reflejan zonas de escaso cambio en la 
intensidad de retrodispersión como es el caso de las 
zonas urbanas (San Pedro, Baradero, Gueleguay) 
que presentan alta intensidad en ambas imágenes. 
En la imagen LANDSAT (Fig. 3) se observa en rojo 
(mayor importancia relativa de la banda 7) las zonas 
afectadas por fuego, denotando el suelo desprovisto 
de vegetación y con presencia de carbón y cenizas. 
Fuera de los límites del delta, donde el fuego no 
alcanzó a extenderse, se observan los cultivos y la 
vegetación natural en tonos más claros. 

La estimación de NBR (Fig. 4.b) permite observar 
zonas oscuras (bajo valor NBR) en una gran 
proporción del delta que denotan zonas dañadas por 
el fuego. Se observa una gran asociación entre la 
imagen de NBR (Fig. 4.b) y la imagen de la diferencia 
de intensidad en el canal HV de ALOS (Fig. 4.a). 
Valores más oscuros de la imagen de NBR se 
corresponden con sectores más oscuros de la 
imagenALOS. 

La comparación visual entre los resultados de la 
clasificación de área quemada a partir de imágenes 
ALOS y LANDSAT (Fig. 5) muestra también cierta 
asociación. La mayor proporción de superficie 
quemada se observa para ambos casos en la zona 
central de la imagen, siguiendo la dirección del río 
Paraná. También se observa en ambas imágenes, 
zonas con baja proporción de áreas quemadas en 
regiones ubicadas al NORESTE y al ESTE. Sin 
embargo, se observan algunas zonas clasificadas por 
ALOS como quemadas, pero no por LANDSAT. 

El análisis cuantitativo de ambas clasificaciones 
(Cuadro 2 — Fig. 6), muestra una asociación entre 
ALOS y LANDSAT superior al 70 % (área clasificada a 
través de ambos sensores como quemada y no 
quemada). Este análisis refleja también mayores 
errores de comisión que de omisión de ALOS, 
mostrando tendencia a sobreestimar el área 
quemada en comparación con la imagen LANDSAT. 
La superficie detectada como quemada por ALOS fue 
del 46.33 %, mientras que la detectada por LANDSAT 
fue del 34.85 %. 
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Figura 4 a) ALOS PALSAR: Diferencia en intensidad HV 2007 vs. 2008. Escala de Grises. Blanco: cambios negativos; 
Negro: cambios positivos (disminución señal HV); Gris: sin cambios. b) LANDSAT TM: indicador NBR: Escala de 
Grises. Blanco: Valores altos; Negro: valores bajos (mayor daño). 


a) Bb) 





j Figura 5. a) Área quemada Landsat TM 6 de Octubre del 2008. 
b) Area quemada Imagen ALOS 2007/2008. Colores: rojo: quemado; gris No Ouemado; Blanco: no clasificado. 


2 1 Edición Especial Vol. 2 


piciemBre 2010. EVALUACIÓN DEL USO DE IMÁGENES DE RADAR EN BANDA L 


REVISTA 


SELPER 


Cuadro 2: Comparación de las clasificaciones de de Área quemada utilizando LANDSAT TM y ALOS PALSAR. 








Resultado comparación Nro. de pixeles 
ACIEATO QUEMADO 3066260 
OMISION 1296507 
COMISION 2976433 
ACIERTO NO OJEMADO 6656772 
Total 14819972 





rlicie [ha rficie 15% 
136720 26.10 
45824 8.75 
105977 20.23 
235277 44 92 
523797 100.00 





Figura 6: Comparación de ALOS vs TM. Colores: Rojo: Acierto quemado; Verde: Omisión; Azul: Comisión; Gris: 
Acierto no quemado. 


Es de destacar que para la clasificación a partir de la 
imagen de RADAR, se consideró una clasificación no 
supervisada en tres categorías (diferencia positiva — 
quemado, diferencia negativa y sin cambios). La 
división que se da entre estas categorías puede no 
necesariamente estar asociada a la identificación de 
daño por fuego. 

La identificación de umbrales que definan estas 
categorías puede ayudar en este sentido. 

La comparación entre focos de calor acumulados 
durante el período considerado, y las clasificaciones 
con imágenes LANDSAT (Cuadro 3) y ALOS (cuadro 
4), mostró cierta asociación, aunque la superficie 
estimada como quemada por Hotspots fue en todos 
los casos menor, como es esperable, dado que sólo 
identifica los focos de calor presentes en el momento 
de la captura de la imagen. El resultado de la 


comparación con ALOS presentó mayor proporción 
de errores que el de la comparación con LANDSAT. 

El cuadro 5 muestra el resultado de las clasificaciones 
según el tipo de vegetación en la zona de estudio. Se 
observa que el resultado de la clasificación varió con 
el tipo de vegetación, mostrando diferencias 
importantes entre ambos sensores en algunos casos 
(praderas de herbáceas acuáticas y bosques) y 
similitudes en otros (juncales). En el caso de praderas 
de herbáceas acuáticas, es posible que la gran 
influencia del agua en estos sistemas pueda influir de 
manera importante en la señal de RADAR y confundir 
la determinación de daño por fuego. Es de destacar 
también que la proporción de área quemada fue muy 
variable entre tipos de vegetación, observándose 
regiones con muy poca área quemada como los 
pastizales de gramíneas. 
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Cuadro 3. Comparación de las clasificaciones de de Área quemada utilizando LANDSAT TM y Hotspots 








Resultado comparación Nro. de pixeles Superficie (ha) Superficie (%) 
ACIERTO QUEMADO 1665124 98552 11,24 
OMISION 3499952 123702 23,64 
COMISION 675456 23674 4.56 
ACIERTO NO QUEMADO 5967673 316953 60.56 
Total 14308715 57333) 100.00 





Cuadro 4. Comparación de las clasificaciones de de Área quemada utilizando ALOS y Hotspots 








Resultado comparación Nro. de pixeles Superficie (ha) Superficie (%) 
ACIERTO QUEMADO 1538936 54392 10.39 
OMISION 5337154 188636 36.04 
COMISION 802279 28156 5.42 
ACIERTO NO QUEMADO 7132143 252078 48.16 
Total 14310512 523463 100.00 





Cuadro 5. Comparación de las clasificaciones de de área quemada utilizando imágenes LANDSAT TM y ALOS 
PALSAR considerando diferentes tipos de vegetación. COR: Pajonales, PRA: Praderas de herbáceas acuáticas 
graminotides y latifoliadas; FOR: bosques; JUN: juncales; GRA: Praderas de gramíneas. 
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El análisis de la clasificación de área quemada con 
ALOS PALSAR utilizando umbrales de cambio en el 
nivel de retrodispersión permitió observar que el % de 
aciertos (menor nivel de errores), se incrementa 
hasta un valor óptimo y luego disminuye (Fig. 7). El 
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umbral óptimo sería cercano a una diferencia de 5 dB, 
observándose para este valor, mayor asociación 
entre el área quemada identificada con ALOS y con 
LANDSAT, y un mayor equilibrio entre la proporción 
de errores de omisión y comisión. 


—— Aciertos 


—$— Quemado 
LANDSAT 


—4— Quemado 
ALOS 


—e— Omisión 
—+— Comisión 


Figura 7. Resultados de la comparación entre clasificaciones LANDSAT y ALOS en función del umbral de 
disminución en el coeficiente de retrodispersión de la señal de RADAR. 


A su vez, considerando cada tipo de vegetación por 
separado se podrían identificar diferentes umbrales 
de cambio en la intensidad de retrodispersión (Fig. 8). 
Los tipos de vegetación con menor nivel de biomasa 
suelen presentar valores de retrodispersión más 
bajos, y por lo tanto, la diferencia entresuelo vegetado 
y suelo desnudo suele ser menor que para el caso de 
tipos de vegetación con alta biomasa (bosques). 

Por esto es esperable que algunos tipos de 
vegetación como pajonales y juncales requieran 
umbrales más bajos para la detección de área como 
quemada que otras de mayor biomasa como los 
bosques. Lo señalado anteriormente en el caso de las 
praderas de herbáceas acuáticas puede generar 
inconvenientes también a la hora de identificar 
umbrales. Con respecto a los pastizales, la muy baja 
proporción de área quemada puede generar ruido 
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también al intentar optimizar un umbral de área 
quemada. 

Como resumen se muestra el resultado de la 
comparación entre imágenes ALOS PALSAR y 
LANDSAT considerando diferentes umbrales de 
cambio en el valor de retrodispersión para cada tipo 
de vegetación (Cuadro 6, Fig. 9). Se observa en este 
caso niveles de acierto mayores al 77 %, y un mayor 
equilibrio entre los errores de omisión y comisión. 
Cabe destacar que aún existen diferencias entre las 
clasificaciones con ALOS y LANDSAT, seguramente 
debido a los efectos confundidos al clasificar con 
ALOS que pueden ser provocados por diferentes 
niveles de humedad del suelo o del canopeo, o debido 
a cambios en la dinámica de la vegetación no 
detectados. 
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Figura 8. Resultados de la comparación entre clasificaciones LANDSAT y ALOS en función del umbral de 
disminución en el coeficiente de retrodispersión de la señal de RADAR para cada tipo de vegetación. COR: 
Pajonales, PRA: Praderas de herbáceas acuáticas graminoides y latifoliadas; FOR: bosques; JUN: juncales; GRA: 
Praderas de gramíneas. 


Cuadro 6. Comparación de las clasificaciones de de área quemada utilizando LANDSAT TM y ALOS PALSAR 





Resultado comparación E a e Superficie (%) 
pixeles (ha) 
ACIERTO QUEMADO 3643608 128780 24.58 
OMISION 1521343 33770 10.26 
COMISION 18/3362 66212 12,64 
ACIERTO NO QUEMADO 77853969 275187 32.52 


total 14824282 323949 100 
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Figura. 9: Comparación de ALOS vs. LANDSAT TM utilizando una clasificación de ALOS con umbral diferencial por 
tipo de vegetación. Colores: Rojo: Acierto quemado; Verde: Omisión; Azul: Comisión; Gris: Acierto no quemado. 


CONCLUSIONES Y RECOMENDACIONES 


En este trabajo se analizó el uso de nuevas 
tecnologías como las imágenes polarimétricas de 
RADAR en banda L para la estimación de daño por 
fuego en la región del delta del río Paraná. Si bien se 
observó gran asociación entre los resultados 
obtenidos a partir de imágenes ópticas en relación 
con las imágenes de RADAR, también se 
encontraron diferencias entre ambas clasificaciones. 
En este trabajo se sugieren umbrales de cambio en el 
valor de retrodispersión de la señal de RADAR para la 
identificación de áreas quemadas que pueden ser la 
base de próximos estudios con RADAR en banda L, 
que hasta el momento han sido escasos. Sería 
recomendable la realización de nuevos estudios 
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considerando series de imágenes multitemporales y 
con menor intervalo de tiempo entre su captura. 
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Análisis del evolución pos incendio de un área con pasturas naturales 
aplicando geotecnologías 


Uboldi, J. A.”; Michalijos, P. A; Barragán, F. * y Mancino, M. * 


" Dpto. de Geografía y Turismo — Universidad Nacional del Sur 


RESUMEN 


El área de estudio comprende un sector del Parque 
Provincial Ernesto Tornquist, ubicado en las Sierras 
Australes, al sud oeste de la Prov. de Buenos Alres. 

El objetivo principal del trabajo fue el de detectar la 
evolución del pastizal natural utilizando imágenes 
satelitales, con motivo de la ocurrencia de un incendio 
intencional el día 20 de enero de 2008, el que se 
extendió hasta el día 24 del mismo mes. El área 
incendiada abarcó 1.600 hectáreas, de un total de 
6.000 has. que posee el Parque, la cobertura que mas 
se vio afectada fue la del pastizal y en mucho menor 
grado la forestal. 

Se utilizaron imágenes obtenidas por el sensor 


Thematic Mapper (226-87), a bordo del satélite 


Landsat V, obtenidas los días 9 de febrero y 12 de 
marzo de 2008, sobre las que se aplicaron diversas 
técnicas digitales: clasificaciones supervisadas, 
índices de vegetación, análisis de componentes 
principales, muestreo por píxeles coligados, además 
de los correspondientes controles de campo. 
Comparando digitalmente y en el campo, la cobertura 
correspondiente al pastizal ubicado en zonas de 
cresta de cerros y en depresiones, se observa que en 
el término de un mes se encuentra en franca 
recuperación, en la imagen del 12 de febrero no se 
aprecian los deterioros ocasionados por el fuego, 
realizando los controles de campo, el daño es 
mínimo. 

La rápida regeneración del pastizal no solo recupera 
el valor paisajístico del área afectada, sino que 
también evita la erosión hídrica en las laderas de las 
sierras y provee de pasturas tiernas a los animales. 


Palabras clave: Pasturas quemadas, Imágenes 
satelitales, Estudio multitemporal, Análisis digital 


ABSTRACT 


The study area includes a sector ofthe Provincial Park 
Ernesto Tornquist, located in the Austral Mountain 
ranges, at the sud west of the Buenos Aires province. 
The primary target of the work is detect the evolution 
of the natural pasture using satellite images, on the 


occasion of the occurrence of an intentional fire the 
20th of January of 2008, the one that extended until 
24th of the same month. The burned area included 
1,600 hectares, of a total of 6,000. The cover that was 
affected was the pasture and in much smaller degree 
the forest one. Images obtained by the sensor 
Thematic Mapper on board the satellite Landsat V 


(226-87) were used, the chosen dates were 9 of 


February and 12 of March of 2008, diverse digital 
techniques were applied: supervised classifications, 
vegetation index, analysis of main components, 
sampling by related pixels, besides the corresponding 
controls of field. Comparing digitally and in the field, 
the cover corresponding to the pasture located in 
zones of hill crest and depressions, is observed that in 
the term of a month itis in frank recovery, in the image 
of the 12 of February are not appraised the 
deteriorations caused by the fire, realised the field 
controls, the damage is minimum. The fast 
regeneration of the pasture not only recovers the 
landscaping value of the affected area, but also it 
avoids the hydric erosion in slopes of the mountain 
ranges and provides with tender pastures to the 
animal. 


Keywords: Burned pastures, Satellite Images, Multi 
temporal study, Digital analysis 


INTRODUCCIÓN 


El fuego es un factor ecológico natural que tiene una 
gran influencia en la evolución de las especies y en la 
dinámica de la vegetación de los pastizales naturales. 
Su ocurrencia forma parte de la estabilidad del 
ecosistema, incidiendo sobre distintas propiedades 
del suelo y sobre la biota que se desarrolla en él. Sin 
embargo el fuego adquiere dimensiones más 
dramáticas cuando afecta económicamente el 
patrimonio del hombre. 

El Parque Provincial Tornquist, de 6.000 has, está 
ubicado en el partido homónimo, sobre la ruta 
provincial N* 76 Km 225, dentro del Sistema de 
Sierras Australes de la provincia de Buenos Aires 


(Fig. 1) 
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Fig. 1: Área de estudio 


Pese a su condición de Parque Provincial, 
agricultores que poseen parcelas cercanas, 
aprovechan las pasturas del Parque para alimentar su 
ganado, en consecuencia la presencia del pastizal 
natural adquiere una importancia económica, 
además de la condición de patrimonio natural para la 
región Durante el período estival (seco y cálido), el 
riesgo de ocurrencia de incendios se multiplica, La 
temperatura media es de 20% C, con máximas 
absolutas de 35% C a 40% C. La vegetación se ve 
beneficiada en su desarrollo por inviernos húmedos y 
fríos en donde la temperatura media es de 8* C, la 


mínima absoluta desciende hasta los -10% C, 


produciéndose heladas y nevadas ocasionales, el 
régimen anual de precipitaciones es de 735 mm, 
ocurriendo la mayor cantidad de las lluvias en 
invierno, mientras que en verano existe déficit de 
agua. 

Estas características crean las condiciones propicias 
para la ocurrencia de incendios, el origen del fuego 
puede ser natural (caída de rayos) o antrópico, siendo 
este último el factor desencadenante de la mayor 
cantidad de siniestros ocurridos en la zona, inclusive 
el caso que nos ocupa en el presente trabajo, ocurrido 
entre los días 20 y 24 de Enero de 2008. 


En el periodo 1999-2006, en el partido de Tornquist, 


se registraron 402 incendios que consumieron 86724 
has, siendo su principal causa la negligencia humana, 
y especificamente dentro del Parque, en los últimos 
12 años, hubo 19 incendios que afectaron 8.000 has. 
(SAyDS). 


OBJETIVOS 
El objetivo principal del trabajo es analizar la 


capacidad de regeneración de las pasturas naturales. 
El objetivo secundario es delimitar con exactitud el 
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área incendiada. 
MATERIALES Y MÉTODOS 


Se busca cumplir los objetivos mencionados a partir 
del análisis digital y visual de imágenes satelitales y la 
elaboración de una cartografía temática específica. 
Las ventanas de trabajo se realizaron sobre las 
imágenes obtenidas por el satélite Landsat V, sensor 
TM, obtenidas los días 9 de Febrero (Fig. 2) y 12 de 
Marzo de 2008 (Fig. 3). 

Sobre las áreas de interés delimitadas en las 
imágenes Landsat se realizó el procesamiento digital 
con el objetivo de realzar el área quemada y analizar 
el desarrollo y evolución de la cobertura vegetal. 

El tratamiento digital implicó en primer lugar la 
corrección geométrica, atmosférica y radiométrica de 
las ventanas ubicadas sobre el área de interés. 

Sobre las ventanas originales se realizó en primer 
lugar una combinación RGB 432 (R:0,7a 1,1pu;G:0,6 
a 0,7 uy B: 0,55 a 0,6 1), sobre ambas fechas se 
tomaron muestras de radiancia sobre píxeles 
coligados (mismo píxel en cada ventana de ambas 
fechas), seleccionando dos sectores, uno de ellos 
son las depresiones, sitio en donde se concentra la 
humedad y existe mayor profundidad de suelo, y el 
otro en los sectores altos de las colinas conformadas 
por los esquistos, en donde existe menor humedad y 
roca aflorante, este procedimiento se realiza para 
detectar el crecimiento de la vegetación, 
predominantemente pasturas, ya que de la 
combinación de bandas realizada es posible construir 
la curva espectral de la vegetación y analizar su 
comportamiento en el espectro visible e infrarrojo 
cercano y en las fechas seleccionadas. 
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Figura 2: Ventana de la imagen del 9 de Febrero 
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Figura 3: Ventana de la imagen del 12 de Marzo 


La presencia parcial de nubes, que cubre un sector 
del área incendiada en la imagen del 9 de Febrero, 
motivó la realización digital de productos derivados 
tales como índices de vegetación (Figs. 4 y 5) y 
análisis de componentes principales, tratando de 
minimizar el efecto nuboso y buscando el realce del 
área quemada para realizar su delimitación con la 
mayor exactitud posible, en consecuencia se 
construyeron sobre ambas imágenes índices de 
vegetación útiles en el análisis de la evolución 
multitemporal de la vegetación, se realizó también un 
análisis de componentes principales, tratando de 
realzar el área quemada. 

Los índices espectrales de vegetación son 
transformaciones que se realizan para caracterizar 
las cubiertas vegetales, reforzando la contribución 
espectral (Pinilla, 1995). Estos índices varían entre 1 


y -1, cuanto mayor sea el resultado, mayor vigor 


vegetal presenta la vegetación, característica que se 
traduce en la imagen con tonos de gris muy claros, en 
contraste con aquellos sectores en los que existe 
escasa presencia o ausencia total de vegetación, en 
donde los tonos de grises van de gris medio a negro, 
en este principio se basan la mayoría de los índices 
de vegetación, combinando las bandas roja e 
infrarroja cercano del espectro (Fórmula 1). 


NDVI = (irc=—r) /( irc + r) 


Fórmula 1: Índice Normalizado de vegetación 
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Figura 4: Índice de vegetación realizado sobre la ventana del 9 de Febrero 
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Figura 5: Índice de vegetación realizado sobre la ventana del 28 de Marzo 


El NDVI posee una alta correlación con el índice de 
superficie foliar (LAI), con la biomasa y con la 
actividad fotosintética (Roque Casanova, J.L., 2000), 
siendo un buen indicador de la salud y actividad del 
vegetal, es hoy en día uno de los elementos clave 
para el estudio de los incendios, ya sean de pasturas 
o forestales, tanto en lo correspondiente a la 
estimación del riesgo como a la determinación del 
área quemada. 

El Análisis de Componentes Principales tiene como 
objetivo resumir un grupo amplio de variables en un 
nuevo conjunto, más pequeño, sin perder una parte 
significativa de la información original. En nuestro 
caso el análisis se realizó sobre seis bandas del 
sensor Thematic Mapper (no se consideró la 
correspondiente al infrarrojo térmico). 


RESULTADOS 


Se realiza la comparación píxel a píxel de los valores 
digitales de radiancia correspondientes a las fechas 9 
de Febrero y 12 de Marzo de 2008, la figura 6 surge 
del promedio de los valores muestreados en las 
depresiones existentes en el área incendiada. 

En ella, la curva correspondiente a los datos del 9 de 
Febrero es el resultado de una notable ausencia de 
vegetación, siendo más representativa de la curva 
espectral del suelo descubierto. Su carácter plano 
ascendente se debe a que en el intervalo 0,52 a 0,60 
B no existe una notable prominencia correspondiente 
a la presencia de clorofila, entre 0,63 a 0,69 y, no 
existe el valle originado en la absorción de energía 
para producir la clorofila que caracteriza en nuestro 
caso la presencia de pasturas, y finalmente, en el 


3 1 Edición Especial Vol. 2 
DICIEMBRE 2010 


intervalo 0,76 a 0,90 y la curva muestra una dirección 
levemente ascendente, comportamiento que indica 
ausencia de vegetación, ya que esta se caracteriza 
por una alta radiancia en el Infrarrojo cercano 
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Observando la curva obtenida con los valores 
correspondientes al 12 de marzo, esta presenta 
claramente las características representativas de la 
presencia de vegetación con clorofila. 


—+— [mag er 12 de 1er 
=- [mugen 2 do fob 


Figura 6: Muestreo en áreas deprimidas 


La figura 7, muestra las curvas construidas a partir de 
muestras tomadas en los sectores elevados (menor 
presencia de vegetación). Las respuestas obtenidas 
son las esperadas tanto para el 9 de febrero como 
para el 12 de marzo, en el caso de la primera fecha, 
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ausencia total de vegetación, mientras que para la 
segunda fecha presencia de vegetación en 
crecimiento (aunque en menor grado que su 
equivalente de la figura 6). 


-—— Tag on 12ds rar 
—— Tag en 9 de feb. 


Figura 7: Muestreo en sitios elevados 
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Figura 8: Resultado del muestreo realizado sobre los índices de vegetación 


Cabe recordar, que las gramíneas pueden 
recuperarse rápidamente después de un incendio, 
debido a que la yema de crecimiento está en la base 
de la planta. Además al observar un pastizal 
solamente se esta viendo entre un 30 y 40% de la 
biomasa vegetal, el resto está enterrado en el suelo, 
conformando las raíces. 

Es por ello la diferencia de valores analizados entre 
las imágenes del 9 de febrero y 12 de marzo. 

Con respecto a la utilización de los índices de 
vegetación en el análisis de la regeneración del 
pastizal, la apreciación se realiza por un lado a partir 
de la visualización de las imágenes obtenidas para 
las dos fechas consideradas (Fig. 4 y 5), en donde los 
colores fríos (verde, celeste, azul) y el negro indican 
ausencia de vegetación, mientras que los colores 
cálidos (amarillo, naranja y rojo) muestran aquellos 
sectores en donde existen diferentes grados de vigor 
vegetal, por otro lado se compararon los valores 
promedio obtenidos de 25 muestras de píxeles 
coligados en ambas fechas, dando lugar al gráfico de 
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la fig. 8, en donde es evidente que los valores 
correspondientes a la imagen del 12 de Marzo, al ser 
mayores, indican una actividad clorofílica incipiente, 
corroborando los aspectos colorimétricos de los 
índices en la zona del incendio, con los trabajos de 
campo realizados el 9 de Febrero y el 15 de Marzo, en 
donde fue posible comprobar una franca 
recuperación del pastizal. 

Analizada la evolución de la vegetación dentro del 
área quemada, resta determinar la extensión de la 
misma, referido a este aspecto y de la realización de 
estudios previos, los mejores resultados fueron 
obtenidos a través del análisis de componentes 
principales, considerando en particular el cuarto 
componente principal, ya que en la imagen resultante 
el área quemada se manifiesta en tonos de grises 
muy claros, casi blancos, mientras que la formación 
nubosa se diluye y aparece en tonos de grises 
oscuros, contrastando notablemente con el área 
quemada. (Fig. 9) 
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Figura 9: Análisis de componentes principales y delimitación del área quemada 
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Las características mencionadas referidas al cuarto 
componente, se deben indudablemente a que este 
sintetiza las características de las longitudes de onda 
correspondientes al infrarrojo, las que al poseer una 
frecuencia mas corta que las radiaciones visibles, 
poseen mayor penetración en la atmósfera, la 
presencia de suelo al descubierto en relación a su 
entorno, hace que este sector se manifieste en tonos 
claros, este análisis es recomendable realizarlo antes 
de efectuar la corrección geométrica de la imagen. El 
área incendiada comprende un total de 1.600 has. 


CONCLUSIONES 


Es evidente la utilidad de las geotecnologías en la 
evaluación y evolución de las áreas afectadas por 
incendios, el análisis digital multitemporal de las 
imágenes permitió un seguimiento del desarrollo de la 
cobertura vegetal afectada, predominantemente 
pastizal (Stipa), análisis reforzado por los controles de 
campo realizados, en donde se corroboró el progreso 
detectado en las imágenes, otra especie en franca 
recuperación (semejante a la del pastizal) y que 
lamentablemente tiene características invasivas es la 
retama (Spartium junceum). Con respecto a la 
cobertura forestal, se vieron afectadas especies del 
genero Pinus, cuya capacidad de regeneración es 
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también alta, pero a largo plazo. 

A menos de dos meses después de ocurrido el 
incendio, en la imagen correspondiente al 12 de 
Marzo de 2008, no se detectan indicios del mismo en 
las imágenes (aunque si en los controles de campo). 
Esta rápida recuperación de la vegetación impidió 
que actúen procesos erosivos, principalmente 
hídricos, sobre la tierra descubierta. 

La elaboración de productos especiales, tales como 
los índices de vegetación y el cuarto componente 
principal, complementó el análisis cuantitativo 
referido a los valores de radiancia de los elementos 
pictóricos (píxel) correspondientes a la vegetación, 
destacando aspectos no apreciables a simple vista en 
las imágenes, tales como la regeneración de la 
cobertura vegetal y el límite del área quemada, 
minimizando a su vez los aspectos atmosféricos 
referidos a presencia de nubes en la imagen 
correspondiente al 9 de Febrero. 

Esta delimitación permitió también evaluar la 
proximidad del área incendiada a las estructuras 
construidas, dato no menor ya que pese al despliegue 
de medios para combatir el fuego, determinadas 
construcciones, inclusive la localidad de Villa 
Ventana, corrieron un riesgo concreto de sufrir los 
efectos del fuego. 


Casanova Roque, J.L. 2000 (ed): Remote Sensing 21 st Century: Economic Environmental Applications. Rotterdam: 


Balkema. 


Pinilla Ruiz, C.1995: Elementos de Teledetección. RA-MA, Madrid, España. Secretaría de Ambiente y Desarrollo 
Sustentable (SAyDS). Dirección de Bosques. Programa Nacional de 


Estadística Forestal. 


3 4 Edición Especial Vol. 2 


DICIEMBRE 2010. MONITOREO DE LOS CAMBIOS EN EL COMPORTAMIENTO... 


REVISTA 


SELPER 


Monitoreo de los cambios en el comportamiento de algunas lagunas 
debido a la oscilación climática utilizando imágenes satelitales en la 
puna jujeña 


Maggi, A.E.”; Navone, S. M. *y Kindgard, F. A. * 


" Facultad de Agronomía, Universidad de Buenos Aires 


RESUMEN 


La Puna argentina es una altiplanicie sobre 3500 
m.s.n.m. de clima árido ubicada en el extremo 
noroeste del país. Corresponde a parte de las 
provincias de Jujuy, Salta y Catamarca. La vegetación 
preponderante es de estepa arbustiva y en zonas 
asociadas a cuerpos o cursos de agua se puede 
encontrar humedales. 

Fueron usadas varias imágenes Landsat (MSS, TM y 
ETM+), sin nubes. Las fechas de adquisición han sido 
escogidas para los meses de Agosto y Septiembre 
coincidiendo con el final de las fases de La Niña/El 
Niño. 

Las imagines Path 232/231 Row 075/076 se 
registraron a partir de imagines ortorectificadas. Para 
el mapeo de las lagunas se utilizó el sistemas de 
coordenadas Universal Transversal Mercator UTM 
195, el mismo se aplico atodas las imágenes. 

Los registros de lluvias habían sido provistos por la 
Estación agro meteorológica del INTA Abra Pampa y 
fueron ordenados considerando la estación de lluvias 
monzónicas desde Julio a Junio. 

En este trabajo se estudiaron las relaciones 
existentes entre las oscilaciones climáticas y el 
comportamiento de la superficie del espejo de agua 
de algunas lagunas de la Puna Jujeña. 

Como resultado del estudio se muestra una relación 
directa entre las precipitaciones y el crecimiento y 
disminución de la laguna de Pozuelos. Ello solo se 
observo durante los eventos Fuertes o de extrema 
magnitud, es decir cuando todas las metodologías 
coinciden en considerar eventos El Niño y La Niña 
(WRCC, CDC, CPC y MEl), 1989 Fuerte La Niña y 
1998 Fuerte El Niño respectivamente. 

La validación del modelo que utiliza a las imágenes 
satelitales para cuantificar las superficies de las 
lagunas como un indicador confiable de las 
oscilaciones climáticas, aún no es lo suficientemente 
consistente y requerirá de una mayor serie de 
imágenes. Las escenas deberían elegirse 
considerando las fases El Niño/La Niña definidas a 
partir de la lista del consenso, por los diferentes 
métodos (WRCC, CDC, CPC y ME!]). 


Palabras claves: ENSO, Puna, Landsat, Laguna 
endorreica. 


ABSTRACT 


The Puna in Argentine is a high plateau on 3500 
MASL with a dry climate located in the extreme 
northwest of the country. lt corresponds to part of the 
provinces of Jujuy, Salta and Catamarca. The 
dominant vegetations are a bush steppe and wetlands 
which are in the moist soils near the water sources. 
Several Landsat images (MSS, TM and ETM+), free 
of clouds were used. The dates of acquisition had 
been chosen from August to September at the end of 
La Niña/El Niño phases. 

The images Path 232/231 Row 075/076 were 
registered using orthorectified images as base. The 
Universal Transverse Mercator Coordinate System 
UTM 19S was used for mapping the lagoons, based 
on WGS84. This system was applied on all the 
available data sets. 

The rainfall records have been provided by the 
agrometeorological service of INTA Abra Pampa and 
they were sorted considering the monsoon station 
from July to June. 

In this work, the relationship between the climatic 
oscillations and the behavior of the water surface in 
some lagoons of the Puna Jujeña is studied. 

As a result of this study, a direct relation between 
precipitations and the increase or decrease in the 
surfaces of the Lagoon of Pozuelos is shown. This is 
only observed during the strong events or in extreme 
magnitude, when in almost all methodologies are 
considered for El Niño and Niña events (WRCC, CDC, 
CPC and MEIl), 1989 Strong La Niña and 1998 Strong 
El Niño respectively. 

The validation of the model that using satellite images 
to quantify the surfaces of the lagoons as a reliable 
indicator of the climate oscillations, is not consistent 
enough and should be required a greater series of 
images. The scenes should be chosen according to 
phases of El Niño/La Niña taken from the consensus 
list defined among the different methods (WRCC, 
CDC, CPC and ME!]). 


Keywords: ENSO, Puna, Landsat, Endorheic 
lagoon. 
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INTRODUCCIÓN 


La Puna argentina es una altiplanicie sobre 3500 
m.s.n.m. de clima árido ubicada en el extremo 
noroeste del país. Corresponde a parte de las 
provincias de Jujuy, Salta y Catamarca. La vegetación 
preponderante es de estepa arbustiva y en zonas 
asociadas a cuerpos o cursos de agua se puede 
encontrar humedales. Tiene una fuerte 
estacionalidad de las lluvias que se presentan 
concentradas durante los meses de verano. Existe 
una marcada amplitud térmica con temperaturas, 
características de zonas áridas. 

No existen grandes ríos que recorran la región, 
siendo las cuencas endorreicas con salares o 
lagunas en los piedemontes, el sistema de drenaje 
mas importante. 

Los suelos son poco desarrollados donde se los 
encuentra y un gran proporción de la superficie esta 
cubierta por roca aún no meteorizada. Limita en las 
alturas con la provincia fitogeográfica Alto andina, 
diferente de la primera en el clima y la vegetación. 

Las precipitaciones de la puna, como se ha 
demostrado en otro ecosistemas áridos sufre una 
gran variabilidad interanual de las precipitaciones 
(Ezcurra, 2006). La disponibilidad de agua pluvial es 
fundamental en numerosos procesos que afectan 


desde los procesos físico-químicos hasta la 


productividad primaria neta de los autótrofos, regula 
las poblaciones de herbívoros y determina la 
disponibilidad de agua para los hombres. 

Se han realizado numerosos estudios a partir de 
documentos históricos así como con técnicas de 
estudio paleoclimáticas que demuestran la existencia 
de oscilaciones del clima y particularmente de las 
precipitaciones en el Altiplano a lo largo de la historia 
del planeta. (Prieto y García Herrera, 2008; Quinn y 
Neal, 1983; Ropelewski y Halpert, 1987). 

Considerando solamente el fenómeno ENSO, es 
importante señalar que no existe una única manera 
para identificar las fases (el Southern Oscillation 
Index (SOl), el Niño 3.4 Index, el Niño 3 Index, o el 
Multivariate ENSO Index) y aún no hay consenso 
sobre cuál es el método más adecuado. Una 
consecuencia de esto es que, aunque están claro qué 
períodos corresponden a Niños fuertes, no ocurre lo 
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mismo con Niños moderados o débiles, 
encontrándose divergencias entre las distintas 
fuentes de información. (Report on the IAl Workshop 
on ENSO and Interannual Climate Variability, 1994; 
Informe Final I|Al ENSO —Argentina, 2004). 

En Sudamérica, los Niños de gran intensidad suelen 
ocasionar inundaciones en el sur de Brasil, Paraguay, 


Uruguay y el norte-nordeste de Argentina, inundacio- 


nes repentinas en las llanuras costeras del sur de 
Ecuador y norte de Perú, sequías severas en el 
noreste de Brasil y sequías en la región del altiplano 
andino. 

El evento de El Niño, como fenómeno integrado 
(hidrográfico y climático), por su aperiodicidad (no se 
produce a intervalos de tiempo determinados), 
interrumpe no solo los regímenes de oleaje e 
hidrográfico, sino también todos los órdenes de la 
vida a través de cambios ambientales en los espacios 
considerados. (Carpio, 1999). 

Estos eventos ocasionan severos daños y pérdidas 


en la sociedad, por ejemplo, el de 1982-1983 dejó 


daños en varios países que alcanzaron los 13 billones 
de dólares y unos 2000 muertos (Glantz, 1984). 

La laguna de Guayatayoc (Figura 1) es un cuerpo 
lacustre semi permanente que recibe las aguas de 
toda la cuenca endógena a la que pertenece. Tiene 
asociados sistemas de humedales asociados a 
diferentes zonas de la laguna. Estos frágiles 
ecosistemas son regularmente sometidos a altas 
cargas ganaderas, que combinadas con ciclos secos 
donde la susceptibilidad de las comunidades 
vegetales aumenta, pueden desencadenar procesos 
de desertificación. Las poblaciones migratorias de 
flamencos dependen fuertemente de la duración y 
expansión del espejo de agua que es sustrato del 
zooplancton del que se alimentan cada ano. A su vez, 
la producción ganadera de camélidos y ovinos se 
basa en el forraje natural que producen las estepas y 
pastizales de la cuenca. Esta producción anual de 
forraje y su distribución en el año esta directamente 
relacionada con las precipitaciones. La carga 
ganadera que soportan apropiadamente las estepas 
en años de lluvias promedio es excesiva para los 
años secos y desencadena la desertificación (Vorano 
y Vargas Gil, 2002). 
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Ubicación de las lagunas de Pozuelos y Guayatayo 
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Figura 1. Ubicación de las lagunas al N Pozuelos al S Guayatayoc en el área de estudio 
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Las poblaciones obtienen el agua de bebida de las 
napas freáticas que se recarga anualmente con los 
aportes pluviales. Frecuentemente al final del período 
seco las poblaciones deben ser auxiliadas con agua 
de otras regiones, además de restringir su consumo y 
la posibilidad de irrigación de los pocos cultivos de 
subsistencia que se encuentran en la cuenca. 

La Laguna de Pozuelos (Figura 1) se sitúa sobre una 
planicie a 3.600 m de altura y constituye uno de los 
humedales más importantes a nivel internacional, 
incluida en la lista de la convención de Ramsar. Esta 
es una laguna permanente pero de salinidad 
fluctuante. En la estación seca, quedan al descubierto 
grandes áreas cenagosas y una rica pradera junto a la 
laguna, que es usada por los lugareños para 
pastoreo, en su mayoría de ovejas y llamas, 
especialmente en el invierno austral. La población de 
flamencos asciende a 25.000 ejemplares, 
conviviendo tres especies: el flamenco austral, la 
Parina grande y la Parina chica, estas últimas nativas 
de la Puna. 

A pesar de la alta variabilidad de las precipitaciones y 
la importancia que tiene en muchos procesos, no hay 
estudios que permitan predecir los efectos de las 
variaciones, anticipar los ciclos y tomar medidas 
preventivas a distintos niveles de decisión. 

Existen algunas evidencias sobre las precipitaciones 
y el fenómeno ENSO en otras regiones andinas 
(Bianchi, 2005; Garreud, 2001) de observarse en 
nuestra región de estudio permitirían predecir con 
anticipación la probabilidad de ciclos húmedos o 
secos. Á su vez los mecanismos que relaciona este 
fenómeno de circulación general de la atmósfera con 
las precipitaciones a escala local son más inciertos y 
han sido propuestos diferentes modelos para la 
región andina. 

Tampoco se conoce que interrelación hidrológica 
existen entre años sucesivos. 

Considerando que existe una serie histórica de 
imágenes satelitales que permiten cubrir los 
principales eventos del fenómeno ENSO en las 
últimas décadas se han realizado algunos intentos 
para estudiar los efectos del cambio climático dentro 
de la zona de la Puna jujeña con imágenes satelitales, 
mediante la observación en los cambios de algunas 
geoformas. (Maggi et als, 2005; Kindgard, Navone, y 
Maggi 2008). 

En un contexto de escasa información meteorológica, 
el uso de sensores remotos resulta una herramienta 
promisoria para analizar los últimos treinta anos de 
precipitaciones y su impacto sobre una variable 
respuesta del sistema hidrológico como es la 
expansión total del espejo de agua de las lagunas de 
Guayatayoc y Pozuelos durante la fase La Niña. 

En este trabajo se estudiaron las relaciones 
existentes entre las oscilaciones climáticas y el 
comportamiento de la superficie del espejo de agua 
de algunas lagunas de la Puna Jujeña. 


MATERIALES Y MÉTODOS 
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El fenómeno ENSO no tiene una única definición. Hay 
diferentes criterios para determinar los años Niño, 
Niña o Neutros. Á veces se llega a un consenso como 
el presentado en la tabla 1, con el agregado de si es 
un fenómeno fuerte o débil. (Null, 2007). 

A partir de la información suministrada por esta tabla 
de consenso se realizó la búsqueda de la mayor parte 
de imágenes Landsat MSS, TM y ETM+ sin nubes y 
con fecha coincidente con la finalización de las fases 
de El Niño y La Niña correspondientes al fenómeno 
ENSO. Tabla 1. Se utilizaron imágenes Landsat MSS, 
TM y ETM + porque existen registros desde 1972, la 
distribución de bandas permite delimitar con facilidad 
cuerpos de agua y la resolución espacial es 
apropiada para la escala 1 : 100.000 de este estudio. 
Se escogieron escenas correspondientes a los 
meses de agosto o septiembre ya que estos 
coinciden con el final de la época seca o antes del 
inicio del ciclo monzónico. Las mismas fueron 
provistas por CONAE, INPE y la USGS. Si bien 
hubiese sido interesante utilizar las imágenes de 
septiembre de 1983 y 1995, al no encontrarse 
disponibles, estas no se pudieron incluir en el análisis. 
Tabla 2. 

Las imágenes Path 232/231 Row 075/076 fueron 
registradas a partir de de imágenes orthorectificadas 
del GLCF UTM 195 WGS84. La Laguna de Pozuelos 
al poseer un cuerpo de agua homogéneo permitió el 
uso del algoritmo SAM para clasificar el área con 
agua de las lagunas, mediante el ajuste de las Rule 
Images se logro identificar los cuerpos de agua que 
fueron exportados como shapefile, se vectorizó la 
clasificación y se calculo superficie cubierta por el 
espejo de agua. Para la delimitación de la superficie 
en la Laguna de Guayatayoc y debido a su mayor 
dispersión del cuerpo de agua, se utilizo la 
combinación falso color de las bandas del infrarrojo 
cercano (banda 4) en el canal rojo, y las dos bandas 
del infrarrojo medio (Bandas 5 y 7) en los canales 
verde y azul respectivamente. 

Estas bandas son seleccionadas por las caracte- 
rísticas espectrales de los cuerpos de agua y zonas 
húmedas (Frazier y Page, 2000). Los mismos 
presentan una baja reflectancia en el rango de 0,76 a 
0,9 micrones (banda 4), y nula reflectancia en los 
rangos superiores a 1,3 micrones, donde se 
encuentran la banda 5 (1,55 a 1,75 micrones) y la 
banda 7 (2,08 a 2,3 micrones). Á su vez el suelo 
desnudo, cobertura preponderante de la cuenca, 
presenta valores crecientes de reflectancia a medida 
que las longitudes de onda son mayores. Para las 
imágenes Landsat MSS se utilizó la combinación 
falso color de las de las bandas del infrarrojo cercano 
(banda 7 y 6) en los canales rojo y verde 
respectivamente y la banda del rojo en el canal azul. 
Esta combinación, permite delimitar el cuerpo de 
agua del suelo desnudo y húmedo. 

Esta delimitación del polígono se realizó a través de la 
digitalización de la interpretación por análisis visual 
de las imágenes listadas en el estudio. La superficie 
de los cuerpos de agua se cuantifico usando una 
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extensión del Arc View Gis (XTOOLS), que calcula la Guayatayoc solo se considero el sector norte, 
superficie a partir de polígonos delineados. exceptuándose del análisis el salar ubicado al sur 
Para cuantificar la superficie en la Laguna de debido a su menor ajuste con la estación lluviosa. 


Tabla 1. Años Niño y Niña, ver columna de Consenso. W = episodios cálidos // C = episodios frios // — intensidad débil // + 


intensidad fuerte. Centro del Clima Regional del Oeste (Western Regional Climate Center - WRCC); CDC- Climatic Data 


Center; CPC Climate Prediction Center CENTRO DE PREDICCIONES CLIMATICAS y Multivariate ENSO Index 
(MEI). Fuente: Jan Null, CCM, 2004. 
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Para la cuantificación de la superficie del espejo de 
agua se utilizaron las siguientes imágenes satelitales 
Landsat TM y MSS correspondientes al mes de 
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septiembre de cada año. Se eligió esta época, por ser 
el límite entre la estación seca y el comienzo 


temprano de las lluvias (Buitrago, 2000). 


Tabla 2. Lista de imágenes satelitales utilizadas en este trabajo. 


POZUELOS 
SENSOR FECHA 
Landsat 3 MSS 4-Ago-1978 
Landsat 5 TM 10-5ep-1985 
Landsat 5 TM 19-5ep-1998 
Landsat 5 TM 21-Ago-1999 
Landsat 5 TM 2¿1-Sep-2001 
Landsat 5 TM 17-5ep-2003 


Las precipitaciones correspondientes a los años 
estudiados fueron obtenidas de los registros de la 
estación meteorológica del INTA Abra pampa. La 
misma se ubica a 22” 48' 10.29”de latitud sur y 65% 49' 
27.74” de longitud oeste a 3454 msnm. Esta estación 
se encuentra en forma equidistante y entre ambas 
lagunas, 44 km al norte de Guayatayoc y 47 km al 
sudeste de Pozuelos. Los valores de precipitación se 
incrementan hacia el norte y el este, es decir en 
dirección a la Laguna de Pozuelos y decrecen en 


GUAYATAYOC 
SENSOR FECHA 
Landsat 3 MSS 9-5ep-1978 
Landsat 5 TM 3-5Sep-1989 
Landsat 5 TM 12-Sep-1998 
Landsat 5 TM 21-Ago-1999 
Landsat 7 ETM + 19-Sep-2001 
Landsat 5 TM 12-Oct-2003 


dirección sur y oeste en dirección a la laguna de 
Guayatayoc. Se obtuvieron y ordenaron las 
precipitaciones considerando la estación monzónica. 
Se reprocesaron para representar la estación lluviosa 
del año, de modo que los valores acumulados para el 


año 2000 corresponden a los meses de julio- 


diciembre del 1999 mas los valores registrados en el 
2000 de enero a Junio. De este modo los valores se 
observan el la tabla 3 y figura 2. 


Tabla 3. Precipitaciones anuales totales por ciclo monzónico de julio a junio, de la estación meteorológica INTA 














Abra Pampa. 

AÑO men Adra man 
1971.72 ara 1985-59 13% 
1972.73 3 1989.90 32) 
1973-74 444 1990-91 426 
1974.75 263 1991.92 220 
197570 305 195972533 ¿5 
1976-77 292 1993-94 313 
1977.78 378 1994.95 3130 
1978.79 395 1995-30 319 
1979-80 236 1996-97 471 
1950-51 ¿59 1997398 120 
1981-52 122 1915-99 156 
193253 178 1999-00 ¿71 
1983.84 653 2000-01 3159 
1984-55 122 2001-02 321 
1935-56 923 2002.03 ¿39% 
1186-87 450 2003-04 342 
1987-68 396 





RESULTADOS Y DISCUSIÓN 


Los valores obtenidos en la medición de superficie del 
espejo de agua en la porción Norte de la laguna de 
Guayatayoc, oscilan entre el máximo 5484 ha 
registrado en 2001 y el mínimo 86.8 ha registrado en 
el año 1998. 

El gráfico 4 muestra la relación entre los valores 
acumulados de precipitaciones entre el 16 de 
septiembre del año anterior y el 15 de septiembre del 


año de la medición, reflejando los mm acumulados 
durante la estación lluviosa del año. 

El promedio de precipitaciones para el período 
estudiado fue de 334.4 mm y el desvío estándar de 
108.8 mm. 

En la tabla 4 se observan las estadísticas de las 
precipitaciones de los períodos Niño y Niña y su 
relación con la superficie de las lagunas analizadas 
para los eventos donde hay imágenes disponibles. 
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Abra Pampa - Lluvias por Ciclos Monzonicos 


mm 


1970 4973 4976 1179 


1142 





1985 


Año 


1088 4004 1004 4907 2066 2003 


Figura 2. Para evaluar el efecto de las precipitaciones sobre las variaciones de superficie de la laguna, se uso 
como variable predictiva las lluvias totales acumuladas, siendo la acumulación de agua pluvial de todo el ciclo 
húmedo del monzón, más allá su distribución. Tabla3 y Figura 2. 


Tabla 4. Niños 1978,1998 y 2003 Niñas 1989,1999 y 2001. 


Pp NIÑO Pp NIÑA POZUELOS GUAYATAYOC 
(mm) (mm) NIÑO (has.) NIÑA (has.] NIÑO (has) NIÑA (has. 
. 264,67 3/1,17 502,08 3538,00 599,0 4070,33 
sb 131,49 23,27 671,69 1434,90 638,38 1314,0% 
cv 49,68 6,27 133,78 15,57 116,43 27,28 


Como se desprende de la tabla 4, se observa que los 
años El Niño los desvíos Standard superan a la media 
lo que marca la gran variabilidad de los datos en los 
años de sequía. Los años La Niña muestran una 
menor variabilidad de los datos tanto en las 
precipitaciones como en la superficie ocupada por las 
lagunas. La superficie de la Laguna de Pozuelos 
medido en los eventos La Niña, mostró el menor 
coeficiente de variación. 

Este comportamiento de la laguna aparentemente 
reflejaría la menor dispersión de los datos debido a la 
menor variabilidad entre las precipitaciones 
registradas durante la fase de oscilación 
correspondiente al período húmedo. 

En la tabla 2 y las figuras 4 y 5 se observa que en 
general para los eventos donde existe coincidencia 
(WRCC, CDC, CPC Y ME!) se observa una diferencia 
100 veces mayor de la superficie ocupada por la 
laguna de Pozuelos durante los eventos La Niña. A 
pesar de las diferencias en las escalas temporales, 
estos resultados registrados a partir de imágenes 


satelitales, coincidirían con los obtenidos en trabajos 
de paleontología de las lagunas en la puna 
observándose variaciones de las mismas debido a las 
oscilaciones climáticas del pasado. (Valero Garces, 
2003). 

Cuando se observa la tabla 1 del consenso de todos 
los métodos de calificación del fenómeno ENSO 
desarrollada por los distintos organismos 
internacionales y donde además su intensidad se 
considera fuerte, las lagunas reflejan la oscilación 
climática modificando su superficie. Figuras 4 y 5 En 
los años que se lo califica como fuerte La Niña, se 
observan coincidencias en el comportamiento de 
ambas lagunas analizadas. En el ciclo monzónico 
1988/89 se registraron valores muy elevados, en 
Pozuelos más de 10000has y en Guayatayoc 
superior a las 2500has. Por otra parte para el ciclo 
1997/1998 calificado como fuerte El Niño las 
superficies cubiertas fueron menos de 100 has en 
ambas lagunas Figuras 3, 4 y 5. 
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Figura 3. Superficie de La laguna de Pozuelos. 





Figura 4. Superficie de La laguna de Guayatayoc. 
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Al analizar los valores de precipitación vemos que los 
valores extremos coinciden con fuertes manifes 
taciones del fenómeno ENSO, Fuerte El Niño 120 mm 
y Fuerte La Niña 398 mm, ello sugeriría que estos 
eventos de mayor magnitud se comportarían de 
forma similar y con una oscilación menor que las 
estudiadas en el paleoclima cuya escala temporal es 
más amplia. Tabla 5 y Figura 5. Sin embargo no existe 
una relación directa entre la longitud de la fase El Niño 
y la disminución en la cantidad de lluvia tal como fue 
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observado por (Valero Garcés et al, 2003). Este 
comportamiento también se observa en el caso del 
año 1989, Fuerte La Niña, que para la Laguna de 
Guayatoc se observo un menor crecimiento de la 
laguna respecto de los otros años con evento La Niña 
mientras que en la Laguna de Pozuelos, cuando se 
compara los eventos la Niña, la mayor superficie 
coincide con el año 1989. Estos resultados muestran 
una vez más la gran variabilidad espacial de las 
precipitaciones. Figura 6. 


Tabla 5. Fuertes eventos de El Niño y La Niña 


Fuerte El Niño 
1972/73 
1982/83 
1991/92 
1997/98 


Fuerte La Niña 
1955/56 
1973/74 
1975/76 
1988/89 


Abra Pampa Precipitación Anual 


E) 150 


OE Niño Mia doña 





Precpración (mm) 


Figura 5. Precipitación acumulada en el ciclo monzónico años Niño y Niña. 


4 3 Edición Especial Vol. 2 


DICIEMBRE 2010. MONITOREO DE LOS CAMBIOS EN EL COMPORTAMIENTO... 


Si bien se observan mayores precipitaciones para el 
ciclo monzónico de Julio a Junio durante los eventos 
La Niña existen algunos casos en los que tanto los 
años neutros, como eventos El Niño (año 1978) 
presentan las mayores precipitaciones acumuladas 
con respecto a los Niños 1998 y 2003. Considerando 
que El Niño tiende a asociarse con bajas 
evaporaciones, (Houston, 2005) quien demostró una 
débil asociación entre la evaporación media anual y el 
fenómeno ENSO en el desierto de Atacama, se 
esperaría el crecimiento de las superficies de las 
lagunas. Esto se observa en el caso de Pozuelos que 
en el año 1978 mostró un espejo de agua, 
probablemente como consecuencia de una menor 
evaporación asociada a una mayor precipitación del 
ciclo monzónico. Figuras 3, 4 y 5. Una de las 
particularidades del fenómeno es que cada evento es 
único, en la forma en que se desarrolla, en su 
duración y en la distribución temporal y espacial, 
especialmente, a escala local de sus impactos. 
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Resultados similares encontraron en el informe 
nacional al comparar los valores de precipitación 
registrados en la región del NOA y el fenómeno ENSO 
no siempre se registran los mayores valores 
únicamente durante los años niña. Por ejemplo en los 
años 1984 y 1986 se pueden considerar neutro o débil 
niña, y débil niño respectivamente se observa un 
inusual crecimiento de la laguna de Pozuelos 
alcanzando un pico máximo de 15383 has. en Junio 
de 1984 valores comparables con el evento fuerte de 
La Niña en medido en junio de 1975. Figura 7. 
Observaciones similares fueron registradas por 
(Aceituno y Garreaud, 1995; Vuille, 1999) en la 
Cordillera de los Andes y el altiplano Boliviano. 
Garreaud y Aceituno, 2001 quienes demuestran que 
la cantidad y el signo de las anomalías de las lluvias 
en el Altiplano es altamente dependiente de la 
localización de la zona de anomalía de los vientos y 
por lo tanto requiere estudios de las anomalías de las 
temperaturas en la troposfera. 
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Figura 6 abcd. Imágenes Landsat para de las superficies Pozuelos y Guayatayoc durante eventos extremos. 
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Los factores que mantiene la variabilidad de la 
repuesta son similares a los que interactúan con el 
cuerpo lacustre total. Cuando se disponga de datos 
de la distribución espacial de los datos de 
precipitaciones, las temperaturas registradas durante 
la estación y la dinámica de los vientos se podrán 
mejorar las interpretaciones de los resultados. Del 
análisis comparado entre los índices utilizados en la 
tabla del consenso y los promedios de precipitaciones 
se pueden extraer algunas tendencias que deberían 
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ser verificadas con más cantidad de eventos. Para la 
Laguna de Guayatayoc no se observa una marcada 
diferencia entre lo sucedido durante los casos 
extremos del fenómeno niño/niña y/ó de mayor 
magnitud respecto de otros años niño/niña, por el 
contrario en la Laguna de Pozuelos se observan 
diferencias de comportamiento de la superficie 
ocupada, cuando se compara entre las fases fuetes y 
débiles (Figuras 3,4 y 6). 


Laguna Pozuelos Junio 


Pp jrarn] Sup [Pas] 








3%) 


0%: 004 
e a y 
Wo launs 





Figura 7. Similares crecimiento de la superficie con agua de la Laguna de Pozuelos en diferentes fases del fenómeno 


CONCLUSIONES 


En consecuencia solamente para la Laguna de 
Pozuelos y durante los eventos considerados fuertes 
o de elevada magnitud, donde existe coincidencia 
entre las distintas metodologías para la 
determinación de un evento Niño o Niña (WRCC, 
CDC, CPC y MEl), se observa una relación directa 
entre las precipitaciones y el crecimiento o 
disminución las superficies de la Laguna. 1989 Fuerte 
Niña y 1998 Fuerte Niño respectivamente. 

La validación del modelo de cuantificación de la 
superficies de las lagunas a partir de imágenes 
satelitales como un indicador fiable de las 
oscilaciones de las precipitaciones de la estación 


ENSO. 


lluviosa anterior, aún no es los suficientemente sólido 
y requerirá una serie de imágenes mayor coincidente 
con las fases Niño/Niña surgidos del consenso entre 
los distintos métodos (WRCC, CDC, CPC y ME!l). 

La laguna de Pozuelos fue más sensible respecto de 
la laguna de Guayatayoc, como indicadora de los 
cambios en la superficie cubierta por agua debido a 
las sequías y precipitaciones extremas. 
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Metodología GIS de identificación de campos con quema ilegal de caña 
en la provincia de Tucumán 
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"Área de Sistemas — Dirección General de Catastro — Poder Ejecutivo de Tucumán 


RESUMEN 


En la provincia de Tucumán la quema ilegal de caña 
ocasiona serios trastornos ambientales a los grandes 
conglomerados urbanos y sus habitantes, así 
también como al suministro de las líneas de alta 
tensión que cruzan los campos incendiados. La 
identificación por parte del Gobierno de los 
responsables como los dueños de los campos 
afectados es dificultosa por diversos motivos, entre 
ellos por la gran cantidad de parcelas rurales 
involucradas y la imposibilidad física de determinar 
cuáles son efectivamente los campos en infracción 
por su difícil accesibilidad y desconocimiento de su 
extensión dentro de la llanura tucumana. 

Se presenta una propuesta de metodología basada 
en la utilización de las herramientas GIS de la 
Dirección General de Catastro de la Provincia para 
identificar a los titulares de las parcelas en infracción. 
La misma consiste en la adquisición de las 
coordenadas geográficas desde el sitio del siniestro 
vía aérea o terrestre para su posterior remisión a 
través de telefonía celular en tiempo real o en forma 
diferida a partir de un relevamiento de múltiples 
zonas. La identificación se efectúa aplicando un 
algoritmo de intersección espacial en la base de datos 
de los puntos remitidos con las ortofotos digitales, las 
imágenes lkonos y el sistema parcelario digitalizado 
georreferenciado. De esta manera se proporciona la 
información recabada a la Secretaría de Medio 
Ambiente para la aplicación de las sanciones 
vigentes. 


Palabras clave: Medio Ambiente — GIS — Incendios — 
Detección 


ABSTRACT 


llegal sugar cane farm field fires are causing several 
environmental problems to urban centers at Tucumán 
province, affecting also the electrical pipelines that 
cross the fields under fire. 


Trying to id the people behind this is becomes a 
challenging task to the government since it is quite 
difficult to properly identify the parcels involved and 
sometimes because on the first hand it is hard to get 
into those fields due to the lack of maintenance of rural 
roads and secondly the real extension of the area 
under fire is unknown. 

This paper addresses a methodology by using GIS 
Tools to properly identify the people that own the 
parcels under infraction by using the Dirección 
General de Catastro web services. 

It consists of the acquisition of GPS coordinates from 
the burning field by flying over or directly from taking a 
measure from the field to later relay the info using the 
cellular grid in realtime or simply in offline mode for a 
later processing. The identification process is done by 
applying GIS spatial algorithms based on the 
surveyed points, using the aid of lkonos images and 
the georreferenced parcel system. In this way 
information related to the parcel such as owner name 
and address is provided to the Environmental Agency 
to apply the fine based on current legislation. 


Keywords: Environment-— GIS — Fires — Detection 
INTRODUCCIÓN 


El tradicional inicio de la temporada de zafra en la 
provincia de Tucumán tiene asociado un impacto 
ambiental severo en grandes centros urbanos con 
repercusión incluso en el suministro eléctrico para 
toda la región del NOA e incluso el sistema 
interconectado nacional. Debido a la quema ilegal de 
grandes extensiones de caña de azúcar en las 
instancias previas a su cultivo enormes columnas de 
humo provocan la masiva precipitación de carbonilla 
en casas y Calles de las ciudades ocasionando 
incluso problemas respiratorios a los habitantes.[1][2] 
Paralelamente, la quema debajo de las líneas de alta 
tensión causa el corte abrupto de energía eléctrica no 
sólo a nivel citadino sino a otras provincias del resto 
del país. 
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Figura 1.1. Imágenes de incendios de campos de caña [1][2] 


Pese a ser una actividad penada por la ley, la baja 
incidencia de las multas sumada a la difícil 
identificación de los dueños de los campos donde 
ocurren los siniestros hacen que esta cultura de 
varias décadas no disminuya en la medida que el 
Gobierno pretende. 

Utilizando tecnologías de sensoramiento remoto en 
combinación con un sistema de información 
geográfica se presenta en este paper una 
metodología de identificación eventual de los dueños 
de los campos basada en la detección del lugar con 
coordenadas GPS para su posterior compulsa con la 
base de datos de la Dirección General de Catastro a 
fin de proporcionar toda la información requerida por 
la Dirección de Medio Ambiente de la Provincia. 

Las funciones principales de la DGC que determinan 
qué datos almacenar y mantener se encuentran 
delineadas por la Ley n* 3.907, que la faculta para 
poder realizar, controlar y mantener actualizado el 
catastro territorial ejerciendo el cumplimiento de las 
acciones y tareas de registro y actualización de los 
inmuebles del territorio provincial. 

Además, el artículo 186 del Código Penal establece 
penas que van de de tres a diez años de prisión a 
quien cause incendios intencionales. Dos años atrás 
se modificaron incluso las leyes vinculadas a la 
quema de caña para reducir de veinte a cinco años el 
plan de erradicación gradual para eliminar 
definitivamente la quema de cañaverales en todo el 
ámbito de la Provincia de Tucumán. 

No obstante toda la legislación vigente, las campañas 
y los controles contra la quema de los cañaverales no 


han sido efectivas en los últimos años ya que sólo en 
2008 se produjeron alrededor de tres incendios por 
día durante mayo y octubre, el período que 
comprende la zafra azucarera, totalizando más de 
800 en ese mismo año, contra los 300 focos 
registrados durante 2005. 


DESCRIPCIÓN DEL SISTEMA 


La Dirección General de Catastro (DGC) de la 
Provincia de Tucumán cuenta con un Sistema de 
Información Geográfica territorial georreferenciado 
que abarca la totalidad del territorio provincial. 

Un sistema de información geográfica contiene 
metadatos referenciados geográficamente acerca de 
geografía, bienes y medioambiente. [3] 

El sistema posee una arquitectura basada en 
tecnología de capas o “layers” que contiene 
información alfanumérica asociada a los polígonos de 
parcelas que representan la unidad catastral de 
identificación a través del número de padrón. 

El circuito de alimentación de información se 
encuentra conformado por un motor Oracle 10g 
Spatial Plus, que almacena tanto la información 
alfanumérica como gráfica. La información 
alfanumérica comprende datos de domicilio de 
parcela, nomenclatura catastral, información de 
planos, antecedentes de dominio y valuación fiscal, la 
cual es mantenida a través de un aplicativo con 
tecnología cliente — servidor desarrollado en Visual 
Basic 6. 
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Figura 2.1. Interfase del sistema informático territorial (SIT) [3] 


La información gráfica es mantenida a través de un 
sistema de presentación de datos de digitales por 
parte de los profesionales en formato de Autocad 
(.dwg), los cuales son procesados y administrados 
con Autocad Map para acceder a la instancia gráfica 
del Oracle a través de la vectorización de sus 
archivos. 

Para publicación de datos se utilizan tres 
arquitecturas: 

2.1. Servidor espacial (Mapguide Server) [4]: se 
encarga de distribuir datos espaciales y 
alfanuméricos a través de un servicio basado en la 
web. Es capaz de enviar y recibir peticiones de 


diferentes tipos de datos. 

2.2. Un formato de archivo de mapeo que puede ser 
embebido en la página web (archivo mwf): contiene la 
descripción de las capas de datos a visualizar con su 
correspondiente configuración. 

2.3. Una aplicación basada en la web en la que los 
mapas pueden ser visualizados por un usuario final a 
través de un navegador de Internet. 

La información geográfica se encuentra georre- 
referenciada utilizando el método Gauss Kruger — 
Posgar 98 -— Faja 3. Con esta referencia todos los 
puntos que son relevados son geográficamente 
referenciados en el mapa. 





Figura 2.2. Arquitectura del sistema Mapguide Server [4] 


49 Edición Especial 0.2 METODOLOGIA GIS DE IDENTIFICACION DE CAMPOS... 


METODOLOGÍA 


La estrategia adoptada por la DGC a efectos de 
contribuir con la problemática ambiental de la quema 
de cañas es la identificación de las parcelas 
catastrales sobre los cuales han sido detectados 
focos de incendio intencional. 

Debido a la naturaleza del diseño de la base de datos 
GIS de la DGC, es posible asociar gráficamente un 
polígono o parcela catastral con su correspondiente 
número de padrón, que es un número de cuenta 
llamado también número de partida en otras 
provincias, y de allí obtener toda la información 
alfanumérica que sea relevante a los fines de la 
detección. Dentro de dicho árbol de información, la 
información del responsable fiscal y titulares de 
dominio son los que fueron seleccionados como 
producto final de la detección. 

A efectos de poder materializar dicha metodología, el 
sistema de información geográfica de la DGC posee 
digitalizadas un total de 340.000 parcelas abarcando 
toda la superficie de la provincia. Las mismas fueron 
producto de una restitución urbana y ortofoto digital 
rural que consistió en la creación de un polígono por 
cada parcela catastral relevada, que posee su 
correspondiente identificación de nomenclatura, 
valuación fiscal, antecedentes de dominio y domicilio 
de parcela. 
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Esta restitución se constituyó en una gran ventaja 
para la implementación de la metodología propuesta. 
Dos tipos de detección han sido abordados, 
dependiendo si éstos son en tiempo real (online) o 
diferido (offline) [5]. Ambos poseen dependencias 
tecnológicas cuya resolución inicial permite el 
abordaje de una metodología innovadora y 
perfectible. 


DETECCIÓN OFFLINE 


El método está dimensionado para ser utilizado por la 
policía ambiental u organismo de control del Medio 
Ambiente en la que se puede relevar por oficio, ya sea 
en forma aérea con un helicóptero o terrestre a través 
de camionetas. 

Esta última opción puede tener limitaciones según 
dónde se ubiquen los campos y dependiendo del 
grado de accesibilidad de los mismos. Un gran 
desafío lo constituye en la Provincia el estado de la 
red caminera secundaria y sobre todo la terciaria en lo 
referente al mantenimiento de las mismas. 
Últimamente la detección por denuncia telefónica se 
incrementó notablemente al verse afectadas 
poblaciones que han sufrido la quema intensa de 
cañaverales a metros de barrios, generalmente del 
Instituto de la Vivienda, que se encuentran próximos a 
extensos campos de caña de azúcar. 


Adquisición 
de los puntos GPS 


Conversión de 
LATILON a X Y 


Creación de shapefile de 
puntos relevados 


intersección espacial de la 
capa Parcelas con la capa 
Puntos GPS 


Consuñta SOL en la base 
alfanumérica de los 
padrones identificados 





Figura 3.1. Diagrama de detección offline [5] 
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Dentro del plan de detección bosquejado se efectuó 
una prueba piloto volando en el helicóptero de la 
gobernación con un receptor Trimble Geoexplorer a 
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fin de tomar puntos GPS de los campos que se 
encontraban bajo incendio o con surcos que 
evidencien una quema masiva.[6][7] 





Figura 3.2. Campo con cañaveral en progreso de inicio de incendio visto desde helicóptero [6] y desde un avión a un 
minuto de haber despegado desde Tucumán [7] 


DETECCIÓN ONLINE 


La necesidad de obtención de datos en un tiempo 
reducido, para la identificación parcelaria y sus 
propiedades (como por ejemplo superficie, titulares 
de dominio y los datos de su ubicación, etc.), es uno 
de los factores principales en el desarrollo de esta 
metodología. 

Si bien es cierto que la tecnología avanza cada día 
más, y momento a momento salen al mercado nuevas 
herramientas, también es acertado afirmar que esa 
tecnología que necesitamos, a veces resulta 
inalcanzable económicamente. Es por ello que la 
solución planteada parte de la premisa de no intentar 
alcanzar a la tecnología en su vertiginoso avance, 
sino utilizar la que tenemos disponible al momento y 
ponerla a trabajar para nosotros. 

Consecuentemente, el medio de comunicación 
elegido para la transmisión de la información 
requerida, fue la telefonía celular, de modo tal, que la 
respuesta a los requerimientos fuera casi en tiempo 
real y por sobre todas las cosas, accesible 
económicamente. 

De esta manera, los agentes intervinientes en algún 
siniestro, como lo es la quema indiscriminada de caña 
de azúcar, puede tener los datos del que será el 
responsable de los hechos, en este caso el titular del 
campo, como así también la superficie del mismo y su 
cercanía con líneas de alta tensión o centros urbanos 
en potencial peligro. 

Originalmente, se pensó en el envío de una palabra 
clave que haga referencia a una acción en los 
servidores de catastro, a un número corto (ej. CLAVE 
al 2020); pero para contratar este servicio, se debían 
pautar un mínimo de consumo y a su vez un costo 
agregado a cada mensaje. Este sistema de 
numeración es fácil y práctico para brindar el servicio 
que ofrecemos desde la Repartición, pero encuadra 
dentro de un perfil comercial al cual no quisimos 


acercarnos, y nuevamente, nos aferramos a la idea 
de utilizar solo lo que teníamos en las manos en ese 
momento. 

Cada compañía celular, hasta hace uno o dos años 
atrás, asociaba el número de un teléfono a una casilla 
de correo, por ejemplo, si un número de móvil era el 
154014256, el servicio de telefonía le asociaba la 


casilla de correo  5438140142560EMPRESA- 


PRESTADORA.com.ar, por lo cual se podía enviar y 
recibir correo electrónico a dicha casilla. Gracias a 
este servicio, el agente puede tomar las coordenadas 
de un determinado lugar, incluso con el GPS 
incorporado a algunos celulares, y luego enviarlos a 
la casilla de correo de la Dirección General de 
Catastro. 

Luego de unos segundos, el servidor de Catastro le 
responde de igual manera al celular, con los datos 
requeridos. El sistema estuvo operativo hasta que 
una de las compañías telefónicas dejo de prestar este 
servicio, no así el resto de la competencia. Esto nos 
limitó en cuanto a que no todos tendrían el servicio a 
su disposición. Es por ello que como opción se 
desarrollo una pequeña interfase Web, a la cual se 
pueden conectar todos los celulares por medio de 
tecnología GPRS, para la obtención inmediata de 
datos. 

El GPRS, General Packet Radio Service, es una 
extensión del Sistema Global para Comunicaciones 
Móviles (Global System for Mobile Communications o 
GSM) para la transmisión de datos no conmutada (o 
por paquetes). 

Existe un servicio similar para los teléfonos móviles 
que del sistema |1S-136. Permite velocidades de 


transferencia de 56 a 114 kbps. GPRS se puede 
utilizar para servicios tales como Wireless Application 
Protocol (WAP) , servicio de mensajes cortos (SMS), 
servicio de mensajería multimedia (MMS), Internet y 
para los servicios de comunicación, como el correo 
electrónico y la World Wide Web (WWW). La 


85 1 Edición Especial Vol. 2 
DICIEMBRE 2010 


transferencia de datos de GPRS se cobra por 
volumen de información transmitida (en kilo o 
megabytes), mientras que la comunicación de datos a 
través de conmutación de circuitos tradicionales se 
factura por minuto de tiempo de conexión, 
independientemente de si el usuario utiliza toda la 
capacidad del canal o está en un estado de 
inactividad. GPRS da mejor rendimiento a la 
conmutación de paquetes de servicios, en 
contraposición a la conmutación de circuitos, donde 
una cierta calidad de servicio (QoS) está garantizada 
durante la conexión. Por este motivo, se considera 
más adecuada la conexión conmutada para servicios 
como la voz que requieren un ancho de banda 
constante durante la transmisión, mientras que los 
servicios de paquetes como GPRS se orientan al 
tráfico de datos. [8] 


DETECCIÓN ONLINEATRAVÉS DE SMS: 


1. El agente captura las coordenadas planas con un 
dispositivo GPS. 


2. Redacta un mensaje de texto en el que consignan 
la clave de operación (3 primeros dígitos) seguido de 
un espacio y luego la variable, que en este caso son 
las coordenadas planas XY Gauss Kruger Faja 3. De 
modo tal que el mensaje queda como sigue: 


XYZ 3633146.71,7078462.67 


3. Se envía el mensaje a la casilla de correo de 
catastro (sms(ddgc-tuc.gov.ar) 


4. En el servidor de catastro corre una rutina llamada 
SMSsit, cada 5 segundos, encargada de recibir el 
correo, identificar la palabra clave (en este caso XYZ) 
y enviar la petición al servidor. 
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5. El servidor corre un procedimiento en el que 
intervienen funciones espaciales, las que logran la 
superposición entre el punto recibido (XY) y la 
restitución parcelaria (base cartográfica). Una vez 
capturado el padrón (identificación de la parcela), el 
procedimiento consulta los datos propios de la 
parcela, entre los cuales se encuentra su superficie, 
como también el nombre y apellido del titular de la 
misma. 


6. SMSsit recoge los resultados obtenidos por el 
servidor y los envía a la casilla de correo del celular. 
Todo ello en menos de 1 minuto y al costo de un 
simple mensaje de texto. 


La opción SMS requirió el desarrollo de una rutina 
automática, que pudiera manipular el envió y 
recepción de correo electrónico, y a su vez la 
interpretación del cuerpo del mensaje de texto. 


DETECCIÓN ONLINE GPRS: 


1. El agente se conecta con su móvil vía GPRS (un 
tanto más caro que un simple mensaje de texto) a la 
URLnhttp://200.45.224.98:84/. 


2. Se presiona sobre la opción “Localizar”, y nos pide 
que ingresemos las coordenadas planas XY. 


3. Al presionar el botón enviar, la pantalla siguiente 
contendrá la información requerida en la misma 
pantalla del celular. 


Esta opción GPRS requirió el desarrollo de una 
pequeña pagina Web adaptada para ser visualizada y 
operada desde un teléfono celular. 





Figura. 3.3. Interfase CATWAP de acceso vía celular [8] 
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RESULTADOS 


Durante el sobrevuelo con helicóptero por el sudeste 
de la Provincia se identificaron más de una docena de 
campos con incendio o rastros de incendios previos. 
Dichos puntos o waypoints fueron posteriormente 
convertidos en coordenadas planas x,y a fin de 
subirlos como capa adicional en el GIS de la DGC. 
Una vez convertidos en formato shapefile, los puntos 
son compulsados con los polígonos de parcelas 
rurales utilizando una función del Arcgis 9.3 llamada 
Selección Por Ubicación. Esta función permite 
identificar elementos basados en su ubicación 
relativa a otros elementos. Este tipo de función es 
conocido como consulta espacial. [9] 

El método de inclusión espacial empleado en este 
caso fue el de Completitud (Completely Contain), en 
la que la geometría de la capa de datos fuente debe 
estar contenida en la geometría de la capa destino 
incluyendo sus límites. La capa de puntos GPS 
relevados fueron intersectados espacialmente con la 
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capa de parcelas rurales a efectos de determinar los 
polígonos de las parcelas resultantes. 

En la modalidad online, a partir del SMS se procedió a 
ejecutar el envío del siguiente mensaje de texto: 


XYZ 3582002.3066,/016178.4410 


La repuesta obtenida del servidor de Catastro fue la 
siguiente: 


Padrón: 73545 
Sup: 191542 m2 
Tit: SALOMON, FORTUNATO 


Cabe destacar que la figura 3.5 muestra un campo de 
cañas incendiado donde puede advertirse el paso de 
un ducto de alta tensión sobre el mismo, con una alta 
probabilidad de interrupción del servicio eléctrico 
cuando las llamas se encuentran debajo de los 
cables.[10] 
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Figura 3.4. Intersección espacial de la base cartográfica (parcelas) con los puntos GPS relevados con y sin imagen 
superpuesta [9] 
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Figura 3.5. Posicionamiento del punto remitido vía celular en la base cartográfica[10] 
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CONCLUSIONES 


El empleo de la tecnología de la red celular permite 
utilizar los servicios GIS de la Dirección General de 
Catastro para una detección inmediata o diferida de 
personas responsables de los campos con focos de 
incendio. 

Por otro lado, la manipulación de puntos espaciales y 
su interrelación con el entorno digital nos permiten 
obtener datos de todo tipo, entre ellos: 


- Distancia entre la zona afectada y una línea de alta 
tensión. 


- Población cercana a la zona de riesgo 

-Autoridades policiales y bomberos mas cercanas 

- Historial. Reincidencia en la practica ilegal 

- Identificación del domicilio del responsable 
RECOMENDACIONES 

Los servicios GIS podrían escalarse a medida de las 
necesidades por parte de los organismos de control 
tales como la Secretaría de Medio Ambiente en lo 
referente a la identificación no sólo del titular de la 
parcela en infracción sino también el responsable 


fiscal asociado a quien le llega la boleta del impuesto 
inmobiliario. 
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La posibilidad de ejecutar operaciones espaciales a 
partir de los puntos recibidos permitiría 
eventualmente generar alertas vía email u otro medio 
de comunicación a organismos tales como Defensa 
Civil, Bomberos o Policía en casos que los incendios 
tengan lugar en proximidad de líneas de alta tensión 
para alertar a las autoridades o incluso a la compañía 
de distribución de electricidad. O bien, alertar a 
comunidades de poblaciones cercanas. 

La reciente adquisición de imágenes Ikonos por parte 
de la Dirección General de Catastro comprendió 
únicamente a los polígonos urbanos masivos de la 
Provincia, totalizando un área de 500 km2 de 
superficie. Dicha compra ha permitido actualizar la 
restitución urbana dentro de la campaña de detección 
de mejoras no declaradas por contribuyentes. 

La calidad de los resultados obtenidos en dicha 
actualización permite proyectar la adquisición de 
polígonos de imaginería rurales que permitan 
clasificar de una mejor manera los campos rurales en 
lo referente a los diversos cultivos que son 
sembrados. De esta manera se haría hincapié en la 
caña de azúcar para poder confirmar los estados 
parcelarios a partir de la superposición de dichas 
imágenes con el sistema parcelario georre- 
ferenciado. 

En la figura 6.1. se puede apreciar la proximidad de 
campos de cañas de azúcar con áreas urbanas que 
se verían perjudicadas ante un incendio como los que 
ocurren diariamente. [11] 





Figura 6.1. Imagen Ikonos con restitución parcelaria donde se aprecian las zonas de costura urbanas y rurales. [11] 
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RESUMEN 


Todavía se plantea con frecuencia que la 
teledetección es una tecnología novedosa en los 
estudios de salinidad del suelo y no es menos cierto, 
pero también ya hay experiencia, trabajos y 
resultados importantes. En este trabajo se exponen 
los elementos fundamentales de la metodología 


RADAR-—FAO para la evaluación, monitoreo y manejo 


de desastres naturales en la agricultura y, su 
aplicación a la salinización como un caso de estudio. 


Palabras clave: salinización, desastre, intensidad, 
impacto, RADAR-FAO 


ABSTRACT 


Today often is said that remote sensing is a new 
technology for soil salinity study and these is true, but 
right now there are important experience, works and 


results. Fundamentals features of RADAR-FAO 


methodology for evaluation, monitoring, management 
and its application on the process of soils salinization 
are explained in this paper. 


Keywords: salinization, disaster, intensity, impact, 
RADAR-FAO 


INTRODUCCIÓN 


EL FENÓMENO DE SALINIZACIÓN DESDE LA 
ÓPTICA FÍSICO-QUÍMICA 


Para que ocurra el proceso de salinización en los 
suelos también llamado proceso Solonchak es 
necesario que existan condiciones para la presencia 
de soluciones saturadas con cationes alcalino y 


alcalino-térreos, y el mecanismo de la interacción 


suelo: clima: hombre puede conllevar al proceso de 
salinización del suelo., la formación de los suelos 
Solonchaks tiene relación no solamente con la 
acumulación de las sales, sino también con el factor 
tiempo para que la presencia de las sales solubles 
ejerza su influencia a través del proceso de formación 
de suelos (PFS), dando lugar a la manifestación de 
las propiedades que caracterizan estos suelos. 

Este proceso se manifiesta principalmente por la 
acumulación y redistribución ascendente de las sales 
solubles. En numerosas regiones del mundo la fuente 
principal de las sales solubles proviene de las rocas 
madres. 

Durante la formación del suelo, las sales solubles se 
disuelven por las aguas de lluvia y son lavadas hacia 
las aguas subterráneas, ríos, lagos y mares. La 
salinidad de las aguas marítimas y oceánicas está 
relacionada con la acumulación de las sales solubles 
que provienen de las tierras de nuestro Planeta. 

Por su origen, la acumulación de sal se divide en 
continental y marítima. La acumulación de sal 
continental está condicionada por la concentración de 
sales que se liberan durante el intemperismo en las 
depresiones de diferente tipo y llanuras sin desague. 
Un ejemplo de ello se muestra en la Figura 1. 





Figura 1. Salinización continental 


SO o. ENFOQUES SOBRE EL IMPACTO DE LA SALINIDAD DEL SUELO 


En condiciones continentales con clima árido y 
semiárido, donde las precipitaciones son menores 
que la evapotranspiración, en el caso de un manto 
subterráneo salino cerca de la superficie, ocurre el 
proceso de salinización natural de las tierras en la 
forma siguiente: 

La capa de agua sobresaturada con sales, asciende 
por elevación capilar, alcanzando la superficie del 
suelo. 


ya 





E TA 
- 
4 dr nn á y o Ñ . 
+ e. 2 


y APS pra - 


REVISTA 


SELPER 


El agua se evapora y las sales solubles se precipitan. 
Este mecanismo se mantiene por el flujo capilar del 
manto cerca de la superficie, estableciéndose en 
forma constante como se muestra en la figura 2. De 
esta manera se acumulan las sales fácilmente 
solubles en el agua. Durante la cristalización las sales 
llenan los poros y espacios vacíos del suelo, 
dispersando las partículas y manteniendo en estado 
seco una consistencia friable del suelo 
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Figura 2. Ascensión del agua salinizada desde el manto freático a la superficie del suelo 


UTILIZACIÓN DE LA TELEDETECCIÓN EN LOS 
ESTUDIOS DE LA SALINIDAD DEL SUELO. 


Los método de teledetección aplicados a los estudios 
de la salinidad pueden clasificarse o dividirse en dos 
grandes grupos ( como en muchas otras aplicaciones 
de la teledetección): los que la información se obtiene 
por apreciación visual de las imágenes o a lo sumo 
mediante un indicador de clase y, los que la 
información se obtiene mediante valores de alguna 
magnitud que representa la imagen o elementos 
ocultos tras los que se observan; tal es el caso de la 
salinidad del suelo, o más específico la conductividad 
eléctrica del suelo cubierto por un cultivo o por la 
vegetación. En la figura se muestra un esquema de 
los principales métodos, sensores e imágenes 
utilizados. 

La salinización representa un incremento de riesgos 
medioambientales, especialmente en áreas regadas. 
Esto reclama por el monitoreo para ayudar a tiempo a 
los tomadores de decisión en el manejo, recuperación 
y rehabilitación de la tierra. El monitoreo de la 
salinidad es grandemente realzado cuando se usan 
datos de teledetección, los cuales sin embargo 


necesitan transformación antes que sean usados 
para la identificación y detección del cambio. La 
detección y la estimación de los cambios entre una 
situación actual y una pasada involucran 
incertidumbres sobre la vecindad, naturaleza y 
magnitud de los cambios, que pueden resolverse 
usando el conocimiento de expertos. 

Los mejores resultados de monitoreo se obtienen 
cuando se integran datos de teledetección con datos 
de campo y de laboratorio. En esta consideración, los 
sistemas de información geográfica ofrecen la 
ventaja de integrar datos de diversa naturaleza en 
términos de escala, tiempo, fuente y estructura. Es el 
reto del investigador identificar los indicadores más 
adecuados de la salinidad para un área particular, tal 
que técnicas apropiadas en tierra y de teledetección 
puedan ser aplicadas para extraer la información de 


una manera exacta y costo-efectiva. 


Una de esas técnicas es la de redes neuronales, las 
cuales permiten ajustar la relación entre las variables 
independientes y dependientes sin ceñirse a una 
función matemática específica y por tanto 
incrementar su exactitud. Esto sólo requiere a cambio 
la reiteración del proceso de aprendizaje hasta lograr 
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el ajuste o precisión requerido, de aquí que se 
escogió para desarrollar esta tesis en combinación 
con los métodos de integración de datos vigentes y el 
uso de imágenes multiespectrales de alta resolución, 
que hasta el presente sólo se había usado como 
clasificación de los suelos afectados por la salinidad. 


EL FENÓMENO DE SALINIZACIÓN COMO 
EVENTO DESASTROSO 


Si nos guiáramos exactamente por la definición de 
“evento” expresada en la metodología RADAR 
(Borgia, A., et. al., 2008), no incluiríamos a la 
salinización del suelo entre éstos. Sin embargo, dado 
que el fenómeno encierra elementos de degradación 
del medio, conducentes a fenómenos más 
impactantes aún para el medio ambiente y la 
agricultura en particular, consideramos de 
importancia realizar su caracterización a partir de los 
elementos considerados por la metodología 
mencionada. 

Para la metodología RADAR, un evento es “un 
fenómeno relativamente corto que perdura en el 
tiempo, de una alta amplitud, que es acompañada de 
una alta degradación de los elementos del medio 
ambiente”. Evidentemente, no es el caso de la 
salinización, que no es un fenómeno súbito. Más bien 
es lento y progresivo. 

Tampoco le es característico una amplitud 
desbordada como pudiera ser el caso de un huracán, 
en el que son evidentes los rangos mínimos y 
máximos de la velocidad de los vientos. Sin embargo, 
la salinización es un fenómeno perdurable en el 
tiempo, causante de la degradación del suelo y otros 
elementos medioambientales, razones más que 
suficientes para considerarlo como un fenómeno 
natural a ser considerado en un análisis como el 
presente. 

De igual manera los autores (Borgia, A., et. al., 2008), 
definen lo que para los efectos de la metodología 
RADAR es un desastre en los términos siguientes: 
“Un desastre es la consecuencia general de un 
evento que ocasiona una disrupción significante de la 
vida normal o al menos de pequeños asentamientos 
humanos”.Visto de esta manera, la salinización del 
suelo, no calificaría como desastre. Lo cierto es que el 
fenómeno que nos ocupa, es causante de un impacto 
significativo directo al medio ambiente, en particular a 
la agricultura, por su alta incidencia en la calidad de 
las producciones agrícolas; e indirectamente, a las 
familias productoras, que dejan de percibir las 
ganancias como resultado de la caída de los bienes y 
servicios derivados de sus producciones. 

Aun con estas limitantes, consideramos un reto, 
analizar la salinización del suelo a partir de los 
cánones de la metodología RADAR, y lo hacemos 
con el fin de profundizar en el modelo conceptual de 
éste y de otros fenómenos que se comportan de 
manera lenta y progresiva y llegar a conclusiones que 
permitan generalizar los procedimientos analizados. 
De acuerdo con la metodología RADAR (Borgia, A., 
et. al., 2008), las características principales de un 
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evento desastroso son cuatro: a) duración; b) 
amplitud; c) magnitud; y d) magnitud local. 
Analicemos este punto de vista. 

a) Duración: “Es el intervalo de tiempo entre la 
situación normal antes y después del fenómeno”. 
Obviamente, para el caso de la salinización es un 
tanto difícil situar los límites de cuándo un suelo 
comienza y termina de salinizarse. 

Es decir, no es tan evidente como en otros fenómenos 
que ocurren en un breve lapso de tiempo. No 
obstante, es nuestro criterio de que por los registros 
históricos, es posible acercarse a una fecha que 
podría ser considerada como la situación normal 
(línea base) del territorio. La fecha después del 
fenómeno, podría tener varios límites, por cuanto, el 
fenómeno es progresivo, lo que implica que su 
amplitud de hecho varía con el tiempo en 
dependencia de la intensidad de las precipitaciones y 
de la actividad que introduzca el hombre ara 
recuperar estos suelos. 

b) Amplitud: “Es la medida cuantitativa de la energía 
máxima por unidad de tiempo asociada con el 
evento”. En el caso que nos ocupa, no hay una 
manifestación exógena de energía como es 
característico de otros fenómenos, que al producirse, 
liberan determinados valores de energía. Por 
ejemplo, un terremoto produce determinados rangos 
de energía, asociados a los fenómenos energéticos 
que ocurren en el interior de la tierra y estos valores 
tienen una amplitud dada en función del tiempo. Para 
el caso que nos ocupa, la amplitud estaría dada en 
cómo varía el nivel de conductividad eléctrica del 
suelo en la unidad de tiempo. 

Cc) Magnitud: “Es una medida cuantitativa del total de 
energía de un evento”. Para el caso de salinización 
del suelo, las magnitudes estarían dadas por los 
valores máximos de amplitud registrados, en este 
caso, los valores de conductividad eléctrica del suelo 
para un área de estudio dado. 

d) Magnitud local: “Es la energía cuantitativa 
localmente asociada con el evento; es una función de 
la duración y la amplitud del evento”. Para el caso de 
la salinización, se manifiesta en las magnitudes de 
conductividad eléctrica del suelo en distintos puntos 
del territorio estudiado y para distintos tiempos. 

Como se mencionó antes, a escala internacional 
existe una carencia general de escalas de intensidad 
para la mayoría de los tipos de amenazas. En 
particular no han sido desarrollados trabajos 
consistentes en este sentido en la agricultura para 
identificar la escala de intensidades apropiadas para 
diferentes tipos de eventos desastrosos. 

De ahí que para una rápida evaluación del impacto, 
sea esencial desarrollar escalas estandarizadas de 
intensidad que puedan ser usadas para evaluar el 
impacto en la agricultura, por cada categoría de 
evento destructivo. En la definición de las escalas de 
intensidad, un número preliminar de reglas generales 
deben ser observadas: 

a) Las escalas deben ser simples y fáciles de medir, 
incluso por personas inexpertas. 

b) La definición de los intervalos entre gradaciones 
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subsecuentes de cada escala, debería expresar 
idealmente una correspondencia directa con la 
gradación de la magnitud local del evento (expresada 
en forma de escala lineal o logarítmica). 

Cc) La escala de intensidad puede tener un límite 
inferior de corte (expresado por el daño mínimo), pero 
no un límite superior. 

d) Los componentes del medio que serán usados 
para definir la intensidad en las diferentes 
gradaciones necesitan registrar los eventos sin que 
estén saturados (totalmente destruidos), ni tampoco 
insensibles (de poca degradación o casi nada). De 
hecho, los mismos componentes, pueden o no ser 
usados para diferentes gradaciones. 

e) En cada gradación debería existir un número 
suficientemente grande de alternativas “registrando” 
componentes del medio, para permitir por 
redundancia y comparación de los resultados 
obtenidos en puntos específicos, temporadas, 
climas, composición de suelos y pendientes. 

f) En la definición de cada escala, deberían ser 
establecidas algunos tipos de relaciones que 
permitan transformar de magnitudes (ppm (partes (de 
sal) por millón (de suelos)) o g de sal/g de suelo) a 
intensidades (dS/ím o NDVI) del evento, y de 
intensidades a porcentaje de pérdidas (ecuación de 
Mass y Hoffman). Estas relaciones (las funciones de 
transferencias) pueden ser complejas, ambiguas o 


indefinidas, y conformarían el KB (Conocimiento- 


Base) indicado. Claramente, con el tiempo y la 
experiencia, estas funciones de transferencias 
pueden llegar a ser mejor definidas y cuantita- 
tivamente más sólidas. 

g) Las observaciones remotas (a distancia) de 
eventos específicos deberían ser vinculadas 
directamente a eventos de intensidad. 

Como resultado del estudio sistemático, de la 
investigación y la síntesis de las experiencias de 
trabajos desarrollados a escala nacional e 
internacional sobre el efecto de la salinidad de los 
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suelos en los componentes del medio agrario, se 
pueden plantear las siguientes consideraciones: 


MATERIALES Y MÉTODOS 


ESTIMACIÓN DE LA SALINIDAD 


Para la aplicación de la metodología RADAR-FAO se 
toma como materiales los trabajos realizados en el 
marco del proyecto El-479 Mejoramiento del manejo 


de la salinidad de los suelos cultivados con caña de 
azúcar usando la Geoinformación, coordinado por la 
Universidad de Gante, Bélgica y con la participación 
de la Universidad Agraria de la Habana, Cuba y la 
Universidad Católica de Leuven, Bélgica, vigente de 
1999 a 2004. Los trabajos consistieron en la 
estimación de la salinidad del suelo cultivado con 
caña de azúcar usando una imagen multiespectral 
IKONOS y técnicas de redes neuronales artificiales. 
El sitio experimental se encuentra en la Unidad 
Básica de Producción Cooperativa (UBPC) “Lázaro 
Romero” de la Empresa Azucarera (EA) Héctor 
Molina, en el municipio San Nicolás de Bari, provincia 
La Habana. La UBPC está situada en los 22” 22” de 
latitud norte y los 81* 56” de longitud oeste (ver Fig. 
3.1); tiene una superficie de 1272 ha, de las cuales 
1320 ha están dedicadas al cultivo de la caña de 
azúcar y el resto a cultivos para el autoconsumo, 
principalmente arroz. El clima es tropical, con vientos 
moderados, variaciones de temperatura promedio 
entre 21,3 *C en enero y 27,7 *C en julio; el período 
lluvioso está comprendido entre mayo y noviembre y 
proporciona casi el 75 % del promedio anual de 
precipitaciones que es de 1270 mm. Durante el año la 
humedad relativa fluctúa entre el 60 % y el 90 %. En la 
tabla 1 se muestra un resumen del objetivo del 
trabajo, los lugares donde se hicieron las mediciones, 
los tipos de mediciones y los materiales usados. 


Tabla 1. Resumen de los trabajos realizados para la estimación de la salinidad. 





Lugares de Las Tipo de 
Experimento y 
medicione", medición 
Ajaste de la CFr 
CE: 
con los indices Mogues 15, 
Comdenadis 
de vegetación 2E, 2IN y 
Hurnediard 
mediante redes, 39%. 
Temperatura 
neuronales 


Miteriales, 


imagen multiespectral de satélite IJKLKONOS 
Receptores, GPS 
Sonda de cuatro electrodos (con termómetro 
incorporade)] 
Cilindro porta-muestra, dbalanza, estuta 
Medios de cómputo 
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De manera muy resumida el trabajo consistió en los 
pasos siguientes: 


Toma de datos y preparación de datos 


1. Medición de la conductividad eléctrica aparente 
mediante una soda de cuatro electrodos en 36 puntos 
de 4 bloques escogidos según criterio de experto y 
considerando un desarrollo estratificado de la caña 
de bloques y que tuvieran atenciones culturales o 
agronómicas uniformes. 


2. Determinación de la humedad del suelo por 
métodos de gravimétricos 

3. Obtención de la conductividad eléctrica del extracto 
de saturación a partir de la conductividad eléctrica 
aparente y la humedad. 


Procesamiento de datos 
1. Entrenamiento de la red neuronal con los datos de 


CEe e NDVIl.mediante el programa RSalinidad y la 
realización de números de iteraciones variado. 
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2. Estimación de la salinidad mediante el control de la 
red neuronal Visualización de los mapas obtenidos 
mediante el SIG ILWIS que se muestran en las figuras 
3,4,5y6. 


Los antecedentes de este trabajo de estimación de la 
salinidad del suelo son los trabajos realizados por 
Pulido, etal, 1995 y 1977, Dwivedi, et, al, 1992, 1996 , 
2001., Lau et. al, 1998, 2003, Rao et al, 1991, 1995 y 
Weigand et al, 1991, 1993, 1994, 1996 quienes 
usaron imágenes Landsat SPOT, mediciones de 
campo de conductividad eléctrica, y herramienta 
matemática de regresión lineal simple y múltiple y 
obtuvieron mapas de clasificación y de estimación de 
la salinidad. 

Otros trabajos sobre salinidad de los suelos aplicando 
diferentes técnicas y herramientas pero dentro de la 
teledetección son los realizados por Evan et al, 1995, 
1998 y 2000, Metternicht, et al, 1998a, 1998b, 1999, 
2001., tales como integración de datos, estimación de 
riesgos de salinización, clasificación difusa, relación 
entre reflectancia e indicadores ambientales, 
estimación de cambios espaciales y temporales de la 
salinidad, etc. 
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Figura 3 Mapa de salinidad del bloque 15 
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SALINIDAD DEL BLOQUE 28 
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Figura 4. salinidad del bloqie 38 


ELEMENTOS TOMADOS DE LA METODOLOGÍA 
RADAR 


La aplicación de la ecuación de Maas y Horfman al 
cartograma digital de salinidad nos daría un mapa (un 
raster) de rendimiento del cultivo. La dificultad 
principal aquí es que no se disponen de valores de 
parámetros de la ecuación de Maas y Hofman en la 
mayoría de las situaciones. 

Por último, en el caso que se tuviera tales parámetros, 
al mapa (raster) de rendimiento se realiza un 
agrupamiento de los píxeles (o celdas raster) según 
criterio de experto o de la experiencia (KB) y se 
obtiene un mapa de intensidad del daño de la 
salinización. 

En este caso, dado la similitud de los valores de CE 


dentro de cada bloque, el mapa de intensidad 
aparecería alo sumo de dos a tres colores. 

Como se dijo al principio del trabajo, la salinización es 
un evento desastre bien diferente a un huracán y a un 
terremoto, porque además de no ser brusco ni tener 
tiempo de inicio y conclusión definidos, afecta pocos 
elementos del paisaje en si y mucho menos 
detectables mediante imágenes aeroespaciales: en 
la práctica se reduce al suelo, la vegetación (natural y 
cultivada) y algunos pocos elementos de la 


infraestructura— 


Las Tablas 2, 3 y 4 definen elementos claves 
importantes relacionados con el área donde tiene 
lugar la evaluación del impacto del evento. 


Tabla 2. Área de impacto y comunidad afectada aplicada a GUANAMÓN 





brea Comunidad slectada 


Región Area cañera de la £A “Héctor Molina”. 


Área UBPC "Lázaro Rornero” 


Parcela Bloques cañeros limitrofes con el mangle 


Comunidad Guaramón y aledañas, inciuye hasta los ressdentes de los alrededores del central, 
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Figura 5. Mapa de salinidad del bloque 29N 
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Figura 6. Mapa de salinidad del bloque 298 


APLICACIÓN DE LA METODOLOGÍA RADAR- 
FAO 


La aplicación de la metodología RADAR-—FAO para la 


evaluación de los daños a la agricultura por la 
salinización de los suelos implica o incluye la 
estimación de la salinidad previamente por uno de los 
métodos descritos en el párrafo anterior. 

Para el caso de la salinidad el modelo físico mediante 
el cual se convierte la magnitud (salinidad medida por 
la conductividad eléctrica) en pérdida del rendimiento 
del cultivo es la ecuación de Maas y Hoffman 1977.; 
cuando el suelo está cultivado por supuesto. La 
escala de intensidad del daño es el resultado de la 
experiencia de muchos años. Le ecuación de Maas y 
Hoffman para el cultivo de la caña de azúcar toma los 
parámetros 


RR = 100* (CEo- CEs)/(CEo-—CE100) 
RR = 100* (4.0-CEs)(4.0- 16.0) 


Dado que se acepta que el suelo con CE de 4 S/m o 
menos la caña de azúcar tiene 100 % de rendimiento, 


mientras que para CE de 16 dS/m o más no se da la 
caña de azúcar. 

Desde el punto de vista computacional esta ecuación 
se aplicó a los mapas de salinidad (raster) en el SIG 
ILWIS, mediante la operación MapCalc con lo cual se 
obtienen los mapas de pérdida de rendimiento de la 
caña de azúcar. 

Finalmente, mediante el SIG ILWIS se reagrupan los 
valores de las celdas raster acorde con la escala de 
intensidad de 6 valores definida como se ha reiterado 
sobre la base de la experiencia se obtiene el mapa de 
intensidad del daño. 


RESULTADOS 


Los resultados que se presentan en esta sección 
están afectados o limitados por la dificultad de 
obtenerlos de las instituciones y oficinas que los 
poseen o en ocasiones no existen tales datos, de aquí 
que hemos considerado en llamar al trabajo un 


acercamiento a la metodología RADAR—FAO. 


En la tabla 3 se presentan los elementos del medio 
ambiente afectado por la salinidad y en la tabla 3 las 


Tabla 3. Ambiente y elementos; medio y componentes 


Medioambiente Elementos naturales y socio-económicos (broad) 


Elemento 


Caña de azúcar 


Milieu (medio) Caña de azúcar, suelo, arroyo. 


Componente 


Bloques cañeros 
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Tabla 4. Valores, porcentaje de pérdidas, degradación y disrupción. 

















Valor La caña de azúcar es la materia prima de la producción del azúcar. 
Porcentaje de perdidas No se ha conseguido el dato. 
radación Debido a la veña salinidad el porcentaje areal es insignificante, 
Disrupción Disminución de la producción de azucar 
Disminución de la producción de azúcar (negativa) 
impacto Disminución de la masa foliar 
Cambio de la caña a cultivos para el autoconsumo?. 
impacto total (Toll Disminución de la masa foliar . 
Daño de un com nte Nose tiene el dato 
Daño de una parcela No se tiene el dato 
La figura 7 muestra un esquema de lo que los RADARPFAO aplicada. 


principales elementos de que consta la metodología 


ANÁLISIS DEL MODELO ANÁLISIS EMPÍRICO 


salinnizacion 


EA “Hector Molina” 







Definición de 
Medio ambiente 
impactado 


evento y 
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Ecuacion de Maas y 
Hoffman 


Pocos datos de CE 
Du. 





por componentes 
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Base del conocrmento 
Funciones de transferencia 


$ rro 
33 
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3 SS z Valores de cada 
O Jo componente de 
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INTEGRACION SOBRE LAS coMpONENTES 
DE LAS PARCELAS Y EL ÁREA 


e PARA OBTENER EL DAÑO ) 


Estimación 
del daño 





Figura 7. Elementos de la metodología RADAR-FAO para el caso de la salinizacion. 
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La figura 8 muestra un diagrama de flujo de los 
parámetros a tener en cuenta para caracterizar el 
medio ambiente impactado, y en la tabla 5 se 
presenta la caracterización del evento “salinización 
del suelo” con una estructura de árbol. 


En la tabla 6 se muestra la Serie tentativa de los 
componentes del medio relativos a la producción 
agrícola para los cuales debería ser determinado los 
valores en cada parcela, incluyendo las sub-clases 
para el caso del evento “salinización del suelo”. 





Figura 8. Parámetros de la evaluación del impacto. 


Tabla 5. Evento del “desastre” salinización en estructura de árbol. 





Caracteristicas 


Niveles 
intrusión salina 


Evento Eventos simples causantes 


desastroso 


Mal drenaje 


Niveles de intensidad de los eventos considerados 1-6 





Componentes del medio 
Medio Tipo de cultivos Caña de azúcar 
afectado y su Tipo de suelo Gley Ferralitico Concresionado 
vulnerabilidad Rango de pendientes 14 % 


Tasa de cobertura del suelo 85 % 


Resultados Nivel de daños 


Mecio 
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Tabla 6: Serie tentativa de los componentes del medio relativos a la producción agrícola para los cuales debería ser 
determinado los valores en cada parcela, incluyendo las sub-clases para el caso del evento “salinización del suelo”. 


identíiwcados (de acuerda al media) 
Sistema de recursos para la 
produeción agricola 


A Sistema de ganado No existe 
VIO y ESTBMaS Le sistemas forestales Mangle silvestre] 


Producción 
Construcciones campestre Centros de acopio 
y s5temas ce apoyo Maquinaria y herramienta 
(Incluyendo organización e 
infraestructura) sisterráas de mercado 


Sistemas de investigaci No exsle 





Sistemas de recursas económico Se desconoce 


CONCLUSIONES evaluar los daños en la agricultura por la salinización 
del suelo mediante la teledetección, propicia unificar 

La aplicación de la metodología RADAR-FAO para criterios para evaluar dichos daños, 
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RESUMEN 

Se propone un modelo de predicción de inicio de la 
sequía en el Perú y la costa norte utilizando la serie de 
tiempo de los promedios mensuales de las imágenes 
NDVI - AVHRR - NOAA, los índices El Niño 


Oscilación Sur (ENOS), y la técnica de regresión 
lineal múltiple. Se obtuvieron la serie de tiempo de los 
promedios mensuales de las anomalías 
estandarizadas de los índices NDVI y ENOS en el 
periodo de julio de 1981 a diciembre 2003. La 
evolución mensual del NDVI en el año 2000 en todo 
Perú muestra que los meses de sequía comienzan en 
agosto y terminan en diciembre. El NDVI promedio de 
todo Perú varía de 0.0 (suelo desnudo) a 0.7 (suelo 
vegetado). La desviación estándar del NDVI| muestra 
que en la costa norte varía menos en comparación 
con la sierra y la selva. 


Palabras clave: NDVI, NOAA, Sequía, Índices 
ENOS. 


ABSTRACT 


A model to predict the onset of draught in Peru and in 
the northern coast is proposed which uses the time 


series of monthly averages from the NDVI-AVHRR- 


NOAA images, the El Niño Southern Oscillation 
(ENOS) indices, and the technique of multiple linear 
regression. The time series of monthly averages of 
the standard anomalies of the NDVI and ENOS 
indices for the period from July 1991 to December 
2003 was obtained. The monthly evolution of NDVI for 
the year 2000 in all Peru shows that the months of 
draught start in August and end in December. The 
mean NDV!I in all Peru varies from 0.0 (naked soil) to 
0.7 (soil with vegetation). The standard deviation of 
NDVI shows that in the northern coast it varies less in 
comparison with the Sierra and Jungle. 


Keywords : NDVI, NOAA, Draught, ENOS indices 


INTRODUCCIÓN 


Los avances en el uso de los datos de satélite para 
vigilar las condiciones de la vegetación de la 
superficie de la Tierra hacen posible investigar la 
evolución espacial y temporal de la sequía desde el 
espacio en una escala global y regional. Desde 1981 
los datos colectados por el sensor AVHRR (Advanced 
Very High Resolution) han sido utilizados para 
generar índices tales como el NDVI. Este índice ha 
sido aplicado exitosamente para clasificar los tipos de 
vegetación (Tucker et al., 1985; Menenti et al., 1993), 
y vigilar las condiciones de crecimiento de la 
vegetación (Kogan, 1990; Prince, 1991). El NDVI ha 
sido usado para hacer el seguimiento del estado 
fenológico del crecimiento de la vegetación a una 
escala regional (Justice et al., 1985), para vigilar la 
evolución de la sequía (Kogan y Sullivan, 1993; Liu y 
Kogan, 1996) y estimar la productividad de las 
cosechas (Prince, 1991). Liu y Ferreira (1991) 
reportaron una buena correlación entre la 
precipitación total mensual y el NDVI| mensual 
acumulado. Di, Rundquist, y Han, (1994), a través de 
una simulación numérica, encontraron una buena 
correlación entre la precipitación y el NDVI diarios. 
Por lo tanto, el NDVI| puede ser utilizado 
apropiadamente para inferir la sequía como un 
resultado de la escasez de precipitación. Los datos de 
precipitación colectados por estaciones de 
observación a menudo poseen una pobre resolución 
espacial, especialmente, en regiones remotas de 
difícil acceso, por consiguiente, el NDVI proporciona 
una valiosa fuente de información para vigilar la 
sequía con una alta resolución espacial (Justiniano y 
Luna, 2005; Kogan, 1997). El objetivo de este trabajo 
es contribuir al desarrollo e implementación de un 
procedimiento de predicción de periodos de sequía 
para Perú y la región Costa Norte de Perú usando las 


imágenes NDVI - AVHRR mensual en el periodo julio 
1981 hasta diciembre 2003, los datos TSM medidos 
in-situ de las Regiones Oceánicas 1, 2, 3, y 4, y los 
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datos lOS (Liu y Negron, 2001). 


DATOS IMÁGENES 
DATOS IMÁGENES NDVI - AVHRR-NOAA 


El conjunto de datos GIMMS (Global Inventory 
Modeling and Mapping Studies) (GLCF, 2008; Tucker 
et al. 2004) es un producto denominado índice de 
vegetación de diferencia normalizada (Normalized 
Difference Vegetation Index, NDVI) disponible 
(acceso libre) por un periodo de 22 años desde 1981 a 
2003. Este conjunto de datos es obtenido de las 
imágenes procedentes del instrumento Advanced 
Very High Resolution Radiometer (AVHRR) a bordo 
de la serie de satélites NOAA 7, 9, 11, 14, y 16. El 
conjunto de datos NDVI han sido corregidos por 
calibración, visión geométrica, aerosoles de los 
volcanes y otros efectos no relacionados al cambio de 
vegetación. El NDVI, en general, es una medida del 
estado hídrico de la vegetación y es usado en 
modelos climáticos para calcular fotosíntesis, el 
intercambio de COZ2 entre la atmósfera y la superficie 
del suelo, evapotranspiración y la absorción y emisor 


“4 : 
m4 





METODOLOGÍA 
ÁREADE ESTUDIO 


El área de estudio comprende la región Costa Norte 
(Figura 1a). Se determinara las anomalías NDVI| de 
esa región y se correlacionará con las anomalías 
estandarizada TSM y lOS en el mismo período. 


ANOMALÍA ESTANDARIZADA DEL NDVI 


El NDVI ha sido calculado utilizando la banda 1 (0.58 
— 0.68 um) y la banda 2 (0.725 — 1.10 um) del sensor 
AVHRR de la serie de satélites NOAA, utilizando la 
siguiente ecuación: 


A. <= ¿) 
NDV] = £_2 
p,+p 
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de energía por la superficie del suelo. Los datos 
GIMMS fueron originalmente generados para 
caracterizar cambios biofísicos definidos en el 
proyecto ISLSCP (International Satellite Land 
Surface Climatology Project) (Wilks, 1995). 


Datos índices ENOS 


Los datos usados en este trabajo denominados 
índices ENOS son las anomalías estandarizadas de 
la temperatura de superficie del mar para diferentes 
regiones El Niño y los Indices de Oscilación Sur (10S). 


Estos datos han sido obtenidos de la NOAA (CPC- 


NOAA). En este estudio se han utilizado los datos de 
la Temperatura de la Superficie del Mar de diferentes 
regiones El Niño. La Figura la muestra la ubicación 
de las regiones El Niño. La Temperatura de la 
Superficie del Mar en las regiones oceánicas 1,2 y 3,4 
son obtenidos del Centro de Predicción del Clima 
(CPC-NOAA). El Índice de Oscilación Sur (1OS) es 
definido como la diferencia de presión del nivel del 


mar estandarizado entre Papeete, Tahití (149%W, 
1725) y Darwin, Australia (131*E, 1295). 





(b] 


Figura (La). Mapa del Perú con las ubicaciones de las áreas en estudio. (1b) Regiones El Niño (CPC-NOAA4). 


donde p1 es el valor de la reflectancia de la banda 1, 
p2 es el valor de la reflectancia de la banda 2. El rango 
de variación del NDVI| es de -1 a +1. La serie de 
tiempo de los datos NDVI usados en este estudio son 
los valores promedios de todos los píxeles de 8km x 8 
km de la áreas de estudio (costa norte y Loreto). Estos 
datos son usados para calcular los datos de la 
anomalía de NDVI. Wilks (1995) sugirió que, debido a 
la variación estacional de los eventos climáticos el 
uso de datos de anomalía estandarizada en lugar de 
valores promedio es más conveniente para investigar 
la diferencia NDVI| anual de la normal. La anomalía 
estandarizada, z, es calculada con la siguiente 
ecuación, 


De X, 


O 
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donde z es la anomalía estandarizada, x es el dato 
NDVI, xp es promedio muestral y o es la desviación 
estándar. 


TEMPERATURA DE LA SUPERFICIE DEL MAR Y 
ANOMALIAS DE INDICES ENOS 


El Índice de Oscilación Sur (10S) es definido como la 
diferencia de presión del nivel del mar estandarizado 
entre las estaciones Papeete, Tahití (149%W, 17%5) y 
Darwin, Australia (131%E, 12%5) (GLCF, 2008). En la 
literatura, existen varias definiciones del índice El 
Niño basados en los datos TSM o Anomalías lOS, por 
ejemplo: (a) correlacionada con los valores negativos 
de OS (Walter, 1923; Stoeckenius, 1981; 
Rasmusson, 1983). (b) correlacionada con la TSM 
(Cane et al. 1986) —El evento El Niño está 
caracterizado por un aumento de la TSM más alto que 
2”C por un periodo de 12 meses que resulta en una 
anomalía de TSM en al menos 1*C por al menos tres 


meses; (c) para la región oceánica 3-4, Trenberth 


(1997) definió que el episodio ENOS ocurre mientras 
que el valor de la anomalía TSM es más alto que 
0.4*C (El Niño) y también mientras que el valor de la 
anomalía es mas bajo que -0.4*C (La Niña). En este 
estudio, utilizamos datos TSM de dos regiones 
oceánicas 3,4 (5% N - 52 S y 1202 W -— 170% W) y 1,2 
(0-10%S y 907 W - 80” W) (Figura 1b) 


RESULTADOS 


MAPAS DEL PERÚ: PROMEDIO Y DESVIACIÓN 
ESTÁNDAR DEL NDVI 


En la Figura 2 se muestra la distribución del NDVI 
promedio (2a) y su correspondiente desviación 
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estándar (2b) de todo Perú para la serie de tiempo 
que corresponde a julio 1981 a diciembre 2003 del 


AVHRR - NOAA para píxeles espaciales de 8 x 8 km. 


En la Figura (2a) se observa que el NDVI varía de 0.0 
a 0.7. En la región geográfica costa, el promedio es 
cerca de cero. La región sierra, esta dividida en dos 


sub—regiones, en la del norte el NDVI (0.4-0.5) es 


mayor que en la parte sur (0.1-0.3). En la región selva 
el NDVI es mayor que en las regiones costa y sierra 
(0.5-0.7). La desviación estándar (DE) del NDVI, 


Figura (2b), muestra que en la costa su DE es muy 
pequeño, cerca de 0.0. En la sierra varia de norte a 
sur desde 0.05 a 0.2. En la selva la DE es mayor y 
varia desde 0.15 a 0.3. Ambas Figuras delimitan muy 
bien las tres regiones geográficas peruanas. En la 
Figura 2a, se presenta los límites de las principales 
zonas ecológicas de Brack (2004) sobrepuestas en el 


mapa del NDVI| para el período jul-81-dic-03. La 


heterogeneidad del patrón del NDVI al norte del 
paralelo 9%S se refleja en la multiplicidad de zonas 
ecológicas tipificadas por Brack. Por otro lado, las 
ecozonas siguen el patrón de los Andes en la parte 
central y sur del país. Es evidente que el mapa del 
NDVI contiene información valiosa para afinar la 
posición de los límites de las zonas ecológicas 
propuestas por Brack, en particular, los límites de la 
Selva Alta y Selva Baja, y las zonas al norte del 
paralelo 9.5% S. La Figura 2b muestra los límites de las 
zonas ecológicas de Brack sobre el mapa de la 
desviación estándar (DE) del NDVI. La desviación 
cero corresponde a los desiertos, éstos están bien 
delimitados por el límite de la zona del Desierto del 
Pacífico. La Puna y la Serranía Esteparia muestran 
baja variabilidad. La más alta variabilidad se da en las 
zonas de selva. 


Pr m Pr A 
A bra L a 


NETAS E 
METE RT 


Fig. 2. (a) NDVI promedio: periodo Julio 1981 a diciembre 2003. (b) Desviación Estándar de la serie temporal de NDVI: 
periodo Julio 1981 a diciembre 2003, superpuestas están los límites de las zonas ecológicas de Brack (Brack, 2004) 
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EVOLUCIÓN MENSUAL DEL NDVI EN EL PERÚ 


Los datos NDVI mensuales del Perú fueron 
procesados para el período jul81 a dicO3. En la Figura 
3, se muestra un ejemplo de las imágenes NDV] 
mensuales de Perú para el período de un año 
(ene-dic 2000). La escala de colores indica los 
valores del NDVI, aumentando desde cero, color 
marrón en la izquierda, hasta más alto que 0.7, color 
azul oscuro en la derecha. Las áreas grandes en rojo 
y naranja, que corresponden a valores NDVI más 
bajos que 0.3, aparecieron en agosto a diciembre, 
indicando que la estación seca está centrada en esos 
5 meses. Por lo tanto la iniciación y magnitud de la 
anomalía NDVI negativa dentro de esos cinco meses 
fue usado para caracterizar el comienzo e intensidad 
de la sequía. 


ANOMALÍAS ESTANDARIZADAS DE NDVI, lOS Y 


TSM 
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Los datos de anomalía estandarizada mensual de los 
índices NDVI, lOS y TSM de Julio 1981 a diciembre 
2003 serán usados para desarrollar un procedimiento 
de predicción de la sequía en región costa norte. En la 
Figura 4a y 4b, se muestran los gráficos de la serie de 
tiempo de las anomalías estandarizada NDV!I de la 
Costa Norte versus la anomalía estandarizada de lOS 
y TSM para la región oceánica 1,2 y 3,4. Según la 
comparación de la evolución dinámica del NDVI de la 
costa norte versus la TSM de la región oceánica Niño 
3,4 (Fig. 4) fue observado que, en general, los valores 
de anomalía NDVI responden positivamente a los 
valores de la anomalía de la TSM de las regiones 
oceánicas 3,4, especialmente durante los eventos El 
Niño de 1982/1983, 1986/1987, 1992/1993, 
1997/1998 y 2002/2003, donde fueron observados 
que los picos no son opuestos y en los eventos La 
Niña de 1984/1985, 1988/1989 y 1999/2000, donde 
fueron claramente observados que los picos son 
opuestos. 


Enero Febrero Marzo bril 


Septiembre 


Mayo | Junio 





Julio il higosto | 


Diciembre 


Figura 3. Evolución temporal del NDVI del Perú para el periodo Enero — Diciembre 2000 
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Figura (4a) Anomalía estandarizada del NDVI para la costa norte de Perú y la anomalía estandarizada de la TSM 
de la región oceánica 3.4 para el período Julio 1981 a diciembre 2003. (4b) Anomalía estandarizada del NDVI de 
la costa norte de Perú y la anomalía IOS (Julio 1981 a diciembre 2003). 
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Tabla 1 — Matriz de correlación para la región costa norte. Los valores en negrita son significativamente 
diferentes de O con un nivel de significación a=0.05 
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En la Tabla 1 se muestran los coeficientes de 
correlación del NDVI y los índices ENOS para la 
región Costa Norte, en el periodo 1981-2003. 
Muchos investigadores han considerado que el Clima 
en la Costa Norte de Perú esta influenciado por varios 
sistemas de circulación atmosférica aumentado por la 
TSM en varias regiones y originadas desde el Océano 
Pacifico. Estos sistemas influenciadotes son a 
menudo cuantificados por los índices ENOS. Por lo 
tanto, en este estudio, los índices ENOS con un valor r 
más alto que 0.2 fueron seleccionados como un factor 
del comienzo de la sequía para la construcción del 
modelo. De acuerdo a la Tabla 1, los índices ENOS, 
incluidos TSM de la región oceánica Niño 3,4 y lOS 
fueron usados para desarrollar los modelos del 
comienzo de la sequía ENOS. Con la finalidad de 
verificar que las series de tiempo de cada variable 
presentan un patrón distinto de variación anual cíclica 
o no, los valores de la autocorrelación (r) y los limites 
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de control de las variables han sido obtenidos a partir 
del análisis de autocorrelación de las series de 
tiempo. En la Tabla 2 se muestra el análisis de 
autocorrelación de las series de tiempo, en 
comparación con los coeficientes de correlación (r) y 
los límites de control (+CL) del NDVI, TSM en las 
regiones oceánicas Niño 1+2, Niño 3, Niño 4 y Niño 
3,4; TSM en el Atlántico Norte (NATL), TSM en el 
Atlántico Sur (SATL), TSM en el Trópico (TROP), 
Indice de Oscilación Sur (10S), TSM en la Trayectoria 
1 (SHIP 1), TSM en Dipolo 1 (NATL — SATL) y dipolo 2 
(SATL, NATL). La Tabla 2 muestra que todas las 
variables seleccionadas tienen valores de r mayores 
que el valor del límite de control (Control Limit - CL) al 
menos con cinco meses de desfase en el tiempo (ver 
filas de Lag 1 — 5). Esto demuestra que todas estas 
variables de serie de tiempo tienen una cierto patrón 
de variación cíclico predecible. 


Tabla 2 — Análisis de autocorrelación por comparación del NDVI y los índices ENOS 
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0.675 
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Co] 0381 [0156 | 0.333 ] 0302 | 0019 | 0310 1 0579 | 0333 
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10 | 0,054 | 0.199 | 0.098 | 0.320 | 0.005 | 0.324 | 0293 | 0.370 


o.77s | 0.745 
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Modelo de Sequía ENOS 


En este trabajo estamos considerando que las 
anomalías negativas del NDVI infieren sequía como 
una función del tiempo y de los índices ENOS. Los 
índices ENOS con un coeficiente de correlación (r) 
mayores que 0.2 han sido seleccionados como 
variables independientes. El proceso de correlación 
fue realizado, primero, usando el mismo mes de NDV!| 
e índice ENOS y luego por dislocación el NDVI con un 
desfase de tiempo de un mes con la finalidad de elegir 
el correcto desfase de tiempo para producir el más 
alto coeficiente de correlación. Siguiendo este 
proceso un conjunto de coeficientes de correlación 
fueron obtenidos mientras las series de tiempo de 
anomalías de los índices ENOS de julio 1981 a 
diciembre 2003 se correlacionaron con la serie de las 
anomalías del NDVI de julio 1981 a diciembre de 
2003. Los resultados muestran que los valores de r 
para NDVI/TSM(Niño1+2), NDVI/TSM(NATL), 


NDVI/TSM(SHIP-1) y NDVI/TSM(TROP) fueron 
0.317,0.316, 0.244 y 0.323 respectivamente. 


Modelo 1. El primer modelo ha sido construido 
usando los datos de la anomalía estandarizada del 
NDVI de la costa norte de Perú como variable 
dependiente y las anomalías estandarizadas de la 
TSM de las regiones Niño1+2, SHIP1 y TROP (Figura 
5a) como variables independientes. Los datos 
usados han sido el NDVI| de cuatro meses (Julio a 
diciembre en el periodo 1981 — 1988). Este modelo 
muestra un coeficiente de regresión (R2) de 0.264 a 
un nivel de significancia de 1%. La ecuación del 
Modelo 1 esta dado por: 
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NDVI = -0.146 - 0.551 x TSM(R-Niño1+2) + 0.596 x 
TSM(R-SHIP1)+0.22x TSM(R-TROP) 


Modelo 2. El segundo modelo ha sido construido 
usando los datos de la anomalía estandarizada del 
NDVI de la Costa Norte como variable dependiente y 
la anomalía estandarizada de los índices de la TSM 
de las regiones Niño1+2, SHIP1, TROP y NATL como 
variables independiente. Los datos usados han sido 
el NDVI de seis meses (Julio a diciembre en el periodo 
1981 — 1988). Este modelo muestra un coeficiente de 
regresión (R2) de 0.35 a un nivel de significancia de 
1%. La Figura 5b muestra el modelo observado y el 
simulado. La ecuación del Modelo 2 esta dado por: 


NDVI = -0.123 - 0.578 x TSM(R-Niño1+2) + 
0.664xTSM(R-SHIP1)- 0.007 x TSM(R-TROP) + 
0.332 x TSM(R-NATL) 


Modelo 3. El modelo 3 ha sido construido usando los 
datos de la anomalía estandarizada del NDVI de la 
Costa Norte como variable dependiente y la anomalía 
estandarizada de los índices de la TSM de las 
regiones Niño1+2, SHIP1 y NATL como variable 
independiente. Los datos usados han sido el NDVI de 
seis meses (Julio a diciembre en el periodo 1981 — 
1988). Este modelo muestra un coeficiente de 
regresión (R2) de 0.35 a un nivel de significancia de 
1%. La Figura 5c muestra el modelo observado y el 
simulado. La ecuación del Modelo 3 esta dado por: 


NDVI = -0.123 - 0.580 x TSM(R-Niño1+2) + 0.665 x 
TSM(R-SHIP1)+0.33x TSM(R-NATL) 
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Figura 5a. Modelo 1 de predicción del inicio de la sequía en la Costa Norte. (3b). Modelo 2 de predicción del inicio 
de la sequía en la Costa Norte. (5c). Modelo 3 de predicción del inicio de la sequía en la Costa Norte. 


CONCLUSIONES 


El NDVI promedio, de todos los píxel de 8 x 8 km2, de 
todo Perú varía de 0.0 a 0.7. En la costa el NDVI 
promedio, es muy bajo: cerca de cero. La sierra, esta 


dividida en dos sub—regiones: la parte norte: el NDVI 


es mayor (cerca de 0.4), la parte sur el NDVI| es menor 
(cerca de 0.1). En la selva peruana el NDVI es mayor 
que en las regiones costa y sierra (de 0.5 a 0.7). La 
desviación estándar (DE) del NDVI muestra que en la 
costa la DE es muy pequeño, cerca de 0.0. En la 
sierra varía de norte a sur desde 0.05 a 0.2. En la 
selva la DE es mayor y varía desde 0.15 a 0.3. La 
evolución mensual del NDVI del año 2000 en todo 


Perú presentó baja variabilidad durante todo el año en 
la costa. En la sierra, el NDVI en el norte varió desde 
0.4 a0.5 y en el sur varió desde 0.1 a 0.3. Finalmente, 
el NDVI de la selva varió entre 0.5-0.7. Los valores 
NDVI simulados concuerdan muy bien con los valores 
observados de 1981 a 1989 para los tres modelos. 
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